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Gasification  is  the  thermochemical  conversion  of  solid  fuel  into  the  gas  which  contains  mainly  hydrogen, 
carbon  monoxide,  carbon  dioxide,  methane  and  nitrogen.  In  gasification,  fluidized  bed  technology  is  widely 
used  due  to  its  various  advantageous  features  which  include  high  heat  transfer,  uniform  and  controllable 
temperature  and  favorable  gas-solid  contacting.  Modeling  and  simulation  of  fluidized  bed  gasification  is 
useful  for  optimizing  the  gasifier  design  and  operation  with  minimal  temporal  and  financial  cost.  The 
present  work  investigates  the  different  modeling  approaches  applied  to  the  fluidized  bed  gasification 
systems.  These  models  are  broadly  classified  as  the  equilibrium  model  and  the  rate  based  or  kinetic  model. 
On  the  other  hand,  depending  on  the  description  of  the  hydrodynamic  of  the  bed,  fluidized  bed  models  may 
also  be  classified  as  the  two-phase  flow  model,  the  Euler-Euler  model  and  the  Euler-Lagrange  model. 
Mathematical  formulation  of  each  of  the  model  mentioned  above  and  their  merits  and  demerits  are 
discussed.  Detail  reviews  of  different  model  used  by  different  researchers  with  major  results  obtained  by 
them  are  presented  while  the  special  focus  is  given  on  Euler-Euler  and  Euler-Lagrange  CFD  models. 
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1.  Introduction  and  objective 

The  rapid  growth  in  industrialization  all  over  the  world  with 
simultaneous  increase  in  population,  air  pollution  has  reached  a 
very  critical  level,  which  threatens  the  public  health,  deteriorates 
the  environment  and  damages  property  and  landscape.  An  alarming 
deterioration  of  the  quality  of  the  life  nullified  the  advantages 
gained  by  a  rise  in  living  standards  due  to  industrial  development. 
The  problem  became  so  serious  that  over  the  last  few  decades  all 
industrialized  and  some  developing  countries  introduced  increas¬ 
ingly  stringent  legislation,  restricting  permissible  levels  of  pollutant 
emission  from  major  combustion  systems  such  as  electricity  gen¬ 
erating  power  stations,  furnaces  and  industrial  plants  as  well  as  by 
automobiles  and  aircrafts.  Therefore,  the  conservation  of  limited 
supply  of  fossil  fuel,  climate  change  and  the  increasing  concern  over 
global  warming  prompted  a  search  for  a  new  and  clean  technology. 
Amongst  the  different  technologies,  one  of  the  most  promising 
future  energy  technologies  is  the  fluidized  bed  gasification. 

Gasification  is  the  thermochemical  conversion  of  solid  fuel  into 
the  fuel  gas  which  contains  mainly  hydrogen,  carbon  monoxide, 
carbon  dioxide,  methane  and  nitrogen.  The  product  gas  from  the 
reactor  also  contains  some  contaminants  like  char  particle,  ash  and 
some  higher  hydrocarbons  or  tar.  A  limited  supply  of  oxygen,  air, 
steam  or  a  combination  of  these  serves  as  gasifying  agent.  The 
gasification  consists  of  four  different  steps  e.g.  drying,  pyrolysis  or 
devolatilization,  combustion  or  oxidation  and  gasification  or 
reduction.  These  four  steps  are  described  below. 

2.1.  Drying 

Drying  occurs  at  about  100-200  °C  when  the  moisture  from  the 
solid  fuel  is  driven  out  and  converted  into  vapor.  The  solid  fuel  in 
this  stage  is  not  decomposed  because  the  temperature  is  not  high 
enough  to  cause  any  chemical  reaction. 

1.2.  Pyrolysis  or  devolatilization 

This  is  a  thermal  decomposition  process  where  the  dried  solid 
fuel  is  decomposed  into  low  to  high  molecular  weight  volatiles 
including  tar  and  solid  charcoal  in  the  absence  of  oxygen.  The 
pyrolysis  or  devolatilization  reactions  are  endothermic  and  thus 
the  heat  needed  for  these  reactions  is  supplied  from  the  combus¬ 
tion  or  oxidation  reactions. 

13.  Combustion  or  oxidation 

The  products  of  the  pyrolysis  or  devolatilization  process  are 
partially  oxidized  by  oxygen  supplied  through  air,  and  then  from 
carbon  monoxide,  carbon  dioxide  and  water  vapor  or  steam. 
As  the  combustion  reactions  are  exothermic  and  other  reactions 
in  gasification  are  endothermic,  the  overall  heat  required  for 
endothermic  reactions  is  supplied  by  this  combustion  or  oxidation 
process. 

1.4.  Gasification  or  reduction 

In  gasification  step  several  reduction  reactions  occur  in  absence 
of  oxygen  because  oxygen  is  consumed  in  the  combustion  reactions. 


These  reduction  reactions  are  mostly  endothermic.  The  final  pro¬ 
ducts  from  these  reactions  are  mainly  gas  mixtures  including 
hydrogen,  carbon  monoxide,  carbon  dioxide  and  methane. 

Most  of  the  time  drying,  pyrolysis  or  devolatilization,  combustion 
or  oxidation  and  gasification  or  reduction  steps  are  not  separated  but 
they  are  overlapped  in  gasification  process.  For  example,  in  case  of 
large  particle,  these  steps  can  take  place  simultaneously.  When  the 
large  particle  is  heated  up,  the  outside  portion  is  dried  and  devola- 
tization  also  starts.  In  the  core,  the  particle  is  still  cooler.  When  the 
center  of  a  large  particle  is  heated  up,  too,  and  drying  and  devolatiza- 
tion  already  started  on  the  outside  of  the  particle,  the  residual  char  is 
likely  to  be  already  gasified. 

In  order  to  investigate  the  gasification  process,  different  types 
of  gasifiers  are  developed  like  updraft  gasifier,  downdraft  gasifier, 
cross  draft  gasifier,  bubbling  fluidized  bed  gasifier,  circulating 
fluidized  bed  gasifier,  twin  fluidized  bed  gasifier,  entrained  flow 
gasifier  etc.  Detailed  descriptions  of  the  gasification  technologies 
are  available  in  the  literature  [1,2].  Amongst  the  different  types  of 
gasification  technologies,  the  fluidized  bed  technology  has  a 
number  of  advantages  which  include  but  not  limited  to  the  high 
heat  transfer,  uniform  and  controllable  temperature,  favorable 
gas-solid  contacting,  higher  efficiency  and  fuel  flexibility. 

In  order  to  analyze  the  process  of  fluidized  bed  gasification, 
several  modeling  approaches  have  been  deployed  and  they  are 
broadly  classified  into  two  groups;  equilibrium  modeling  and  rate 
base  or  kinetic  modeling.  Equilibrium  modeling  is  independent  of 
the  type  of  gasifier  because  it  does  not  consider  the  hydrodynamic 
of  the  bed.  Depending  upon  the  process  of  calculating  the  product 
gas  composition,  the  equilibrium  model  may  be  classified  as  the 
stoichiometric  equilibrium  model  or  non-stoichiometric  equili¬ 
brium  model.  Whereas,  the  kinetic  model  generally  consider  the 
hydrodynamics  of  the  bed  coupling  with  the  reaction  kinetics. 
Based  on  the  hydrodynamic  modeling,  the  fluidized  bed  models 
can  also  be  classified  as  two-phase  flow  model,  Euler-Euler  model 
and  Euler-Lagrange  model.  In  hydrodynamic  modeling,  the  most 
established  model  is  the  two-phase  flow  model.  Very  recently,  the 
computational  fluid  dynamic  (CFD)  modeling  of  fluidized  bed 
gasification  based  on  Euler-Euler  approach  and  Euler-Lagrange 
approach  are  attempted  by  researchers  due  to  the  increasing 
computational  power  of  the  modern  computers.  But,  application 
of  CFD  model  to  study  the  fluidized  bed  gasification  process  is  in 
the  developing  stage  and  more  studies  are  needed  [3,4].  In  the 
present  work,  a  detailed  review  of  different  fluidized  bed  gasifica¬ 
tion  models  published  in  the  literature  is  presented.  Mathematical 
equations  governing  fluid  and  solid  flow,  heat  and  mass  transfer 
and  chemical  reactions  for  each  model  are  presented.  Advantages 
and  disadvantages  of  different  modeling  approaches  and  major 
results  obtained  are  discussed.  The  special  attention  has  been 
given  to  the  recently  published  Euler-Euler  and  Euler-Lagrange 
CFD  models. 


2.  Fluidized  bed  gasifiers 

There  are  different  types  of  fluidized  bed  gasifiers  reported  in 
the  literature.  Amongst  them,  a  detailed  description  of  the  bubbling 
fluidized  bed  gasifier,  the  circulating  fluidized  bed  gasifier  and  the 
twin  fluidized  bed  gasifier  are  presented  here.  Fig.  1  shows  the 
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Fig.  1.  Schematic  diagram  of  (A)  bubbling  fluidized  bed  gasifier  and  (B)  circulating  fluidized  bed  gasifier. 


schematic  diagram  of  bubbling,  circulating  and  twin  fluidized  bed 
gasifiers. 

2.1.  Bubbling  fluidized  bed  gasifier 

In  a  typical  gasifier,  when  the  gas  velocity  is  increased, 
a  situation  is  reached  when  the  particles  are  just  suspended  in 
the  upward  flowing  gas.  At  this  situation  the  frictional  force 
between  a  particle  and  fluid  counterbalance  the  weight  of  particle, 
the  vertical  component  of  the  compressive  force  between  adjacent 
particles  disappears  and  the  pressure  drop  through  any  section  of 
the  bed  equals  the  weight  of  fluid  and  particles  in  that  section. 
At  this  instance,  the  bed  is  considered  to  be  at  minimum  fluidiza¬ 
tion  condition.  With  an  increase  in  velocity  beyond  the  minimum 
fluidization  velocity,  large  instabilities  with  bubbling  and  channel¬ 
ing  of  the  gas  are  observed.  At  higher  flow  rates  agitation  becomes 
more  violent  and  the  movement  of  solids  becomes  more  vigorous. 
However,  the  bed  does  not  expand  much  beyond  its  volume  at 
minimum  fluidization.  Such  a  bed  is  called  bubbling  fluidized  bed 
(Fig.  1A).  In  a  bubbling  fluidized  bed,  gas  moves  through  the  bed  in 
void  and  in  bubbles  with  higher  velocity.  Some  of  the  particles  are 
entrained  into  the  freeboard  along  with  these  fast  moving  bubbles 
and  some  fine  particles  are  transported  with  the  product  gas  and 
leave  the  reactor  at  the  top.  But  most  of  the  entrained  particles  fall 
back  and  can  be  continuously  removed  from  the  fluidized  bed 
with  the  remaining  ash  at  the  bottom.  The  bubbling  fluidized  bed 
gasifiers  have  a  number  of  advantages  over  non-fluidized  bed 
gasifiers  which  include  fuel  flexibility,  uniform  temperature, 
higher  efficiency,  lower  capital  and  maintenance  cost  etc. 

2.2.  Circulating  fluidized  bed  gasifier 

When  the  gas  velocity  is  increased  further,  beyond  the  bub¬ 
bling  fluidized  bed  regime,  the  solids  will  be  distributed  across  the 
whole  riser  height  and  entrained  by  the  gas  at  the  top  of  the 
gasifier.  Particles  are  separated  from  the  gas  in  a  cyclone  and  are 
returned  to  the  fluid  bed  near  the  bottom.  Then  it  becomes  a 


circulating  fluidized  bed  (Fig.  IB).  The  advantage  of  circulating 
fluidized  beds  in  mainly  due  to  the  longer  overall  residence  time. 
The  solids  circulate  in  the  outer  loop,  going  up  in  the  riser,  leaving 
at  the  top,  going  down  in  the  return  leg  and  entering  the  riser 
again  at  the  bottom.  But  there  is  also  internal  circulation  of  the 
solids,  which  fall  back  from  the  higher  region  of  the  riser  and 
move  downwards  at  the  riser  wall.  Circulating  fluidized  bed 
gasifiers  are  normally  used  for  large  applications.  It  has  enhanced 
flexibility  over  bubbling  fluidized  bed  gasifier  for  firing  multi-fuels 
with  high  moisture  content  and  significantly  higher  efficiency. 
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2.3.  Twin-fluidized  bed  gasifier 

Another  form  of  fluidized  bed  gasifier  which  is  also  reported  in 
the  literature  is  the  twin  fluidized  bed  gasifier  as  shown  in  Fig.  2. 
In  one  fluidized  bed  (generally  bubbling)  the  gasifying  agent 
(mostly  steam)  is  brought  into  contact  with  the  fuel  fed  to  this 
fluidized  bed.  There  is  continuous  discharge  of  bed  material  with 
unreacted  char  particles  from  this  first  fluidized  bed,  which  is  then 
fed  to  a  second  fluidized  bed  (generally  circulating)  operated  with 
air  and  used  as  a  combustor  to  burn  the  remaining  char  and  to 
heat  the  fluid  bed  particles.  Then  the  hot  bed  material  from  the 
second  fluidized  bed  is  circulated  to  the  first  fluidized  bed  to 
supply  the  heat  for  the  endothermic  steam  gasification  reactions. 
This  technology  is  especially  interesting  for  biomass  gasification 
due  to  the  higher  volatiles  content  in  the  biomass  compared  to 
coal.  Here  only  the  volatiles  are  used  for  synthesis  of  gas  produc¬ 
tion.  The  fixed  carbon  content  is  burned  to  supply  the  heat. 
Therefore,  the  slow  gasification  reactions  of  fixed  carbon  with 
steam  or  carbon  dioxide  are  avoided.  Another  advantage  is  that  if 
steam  is  used  in  the  first  fluidized  bed,  the  producer  gas  is  not 
diluted  by  nitrogen  from  the  air. 


3.  Reaction  mechanism  and  kinetics 

Mathematical  modeling  of  fluidized  bed  systems  can  broadly 
be  classified  into  two  groups;  thermodynamic  equilibrium  model 
and  rate  base  or  kinetic  model  as  described  in  details  below. 


It  requires  a  clearly  defined  reaction  mechanism  that  incorporates 
all  chemical  reactions  and  species  involved.  It  uses  the  elemental 
balance  and  the  equilibrium  constants  of  selected  reactions. 
A  stoichiometric  equilibrium  model  may  also  use  Gibbs  free  energy 
data  to  determine  the  equilibrium  constants  of  a  propose  set  of 
reactions. 

A  brief  description  of  stoichiometric  equilibrium  model  based 
on  the  work  of  Loha  et  al.  [6,7]  is  given  below.  In  stoichiometric 
equilibrium  model,  the  gasification  process  is  represented  by  a 
global  gasification  reaction  as  follows: 

CH*Ov  +  wH20  +  mH20  =  yH2H2  +  ZcoCO +  /c02C02  +  j'/qofhO  +  j'CHyQ-G 

(1) 

where  CHxOy  represents  the  solid  fuel  which  contains  carbon, 
hydrogen  and  oxygen.  Other  elements  are  neglected.  Here,  x  and  y 
represent  numbers  of  atoms  of  hydrogen  and  oxygen  per  number 
of  atom  of  carbon  in  the  solid  fuel  and  can  be  obtained  from  the 
ultimate  analysis  of  the  solid  fuel,  w  is  the  mole  of  moisture  per 
mole  of  dry  ash  free  solid  fuel  and  m  is  the  mole  of  steam  added 
per  mole  of  dry  ash  free  biomass.  On  the  left  hand  side  all  are 
inputs.  On  the  right  hand  side,  jh2 ,  rco.  7co2>  Yh2o  an d  Ych4  are  the 
numbers  of  mole  of  gas  species  on  the  product  gas  and  which  are 
obtained  by  solving  material  balance  and  equilibrium  constant 
equations  of  selected  reactions.  The  material  balance  equations  of 
carbon,  hydrogen  and  oxygen  are  given  below. 

Carbon  balance 

1  =  Yco  +  Yco2  +  Ych 4  (2) 


3.1.  Equilibrium  model 

The  parametric  study  and  the  thermodynamic  analysis  of  the 
gasification  process  by  using  equilibrium  model  is  a  popular 
method  because  it  provides  a  useful  design  aid  in  evaluating  the 
possible  limiting  behavior  of  a  complex  reacting  system  and  also  it 
is  computationally  inexpensive.  Though,  it  is  independent  of 
gasifier  design,  still  it  can  provide  a  guideline  for  process  design, 
evaluation  and  improvement.  It  can  also  be  used  to  study  the 
influence  of  most  important  process  parameters.  For  this  reason 
many  authors  used  the  thermodynamic  equilibrium  mode  to  study 
the  fluidized  bed  gasification  process. 

Equilibrium  models  are  generally  developed  based  on  follow¬ 
ing  assumptions; 

•  The  process  occurs  at  steady  state. 

•  The  gasifier  is  isothermal  and  at  equilibrium  condition. 

•  Reaction  rates  are  fast  enough  and  residence  time  is  sufficient 

to  reach  equilibrium  condition. 

•  Gases  except  H2,  CO,  C02  and  CH4,  N2  and  H20  are  negligible. 

•  Char  contains  only  solid  carbon. 

•  Ash  residue  is  negligible. 

•  The  product  gas  is  at  the  gasifier  temperature. 

•  All  the  gases  obey  the  ideal  gas  law. 

•  Potential  and  kinetic  energies  are  negligible. 

There  are  two  general  approaches  for  equilibrium  modeling  viz. 
stoichiometric  equilibrium  model  and  non-stoichiometric  equilibrium 
model.  However,  both  approaches  are  conceptually  similar  [5],  they 
differs.  The  detailed  description  of  stoichiometric  equilibrium  model 
and  non-stoichiometric  equilibrium  model  is  presented  below. 

3.1.1.  Stoichiometric  equilibrium  model 

The  stoichiometric  equilibrium  model  is  based  on  selecting 
those  species  which  are  present  in  the  largest  amounts,  i.e.  those 
species  having  the  lowest  value  of  Gibbs  free  energy  of  formation. 


Hydrogen  balance 
x+2w+2m  =  2  yHl  +  2yH20  +  4/CHt 


(3) 


Oxygen  balance 

y+w+m  =  Yco +2yC02  +/h2o  (4) 

Now,  it  is  required  to  select  the  major  chemical  reactions  for 
calculating  the  product  gas  composition.  For  the  present  calcula¬ 
tion  following  reactions  are  selected. 

Water-gas  shift  reaction  C0  +  H20  =  C02+H2  (5) 


Methane  reaction  :  C+2H2  =  CH4  (6) 

All  gases  are  assumed  to  be  ideal  and  all  reactions  form  at 
atmospheric  pressure.  Therefore,  the  equilibrium  constants  of  the 
above  two  reactions,  which  are  functions  of  temperature,  are 
given  below. 

The  equilibrium  constant  of  water-gas  shift  reaction  is 
K  =  Pco2Ph2  =  xC02Xh2  7 

PcoPh2o  xcoxh2o 

where  x,  is  the  mole  fraction  of  the  gas  species  iin  the  gas  mixture, 
xi  —  Yi/v total  tmd  here 


Y  total  —  Yh2  +YC0  +  YC02  +  Yh20  +  YcH4- 

The  equilibrium  constant  for  methane  formation  reaction  is 

K  __  pch4  _  *ch4 
2  (Ph2)2  (Xh2)2 


(8) 


Equilibrium  constant  K  is  a  function  of  temperature  and  can 
directly  be  obtained  from  JANAF  thermochemical  tables  [8]  or 
can  be  calculated  from  the  Gibbs  function  as  given  below. 


In  K  = 


A  G° 

Iff 


(9) 


AG£  —  2  /i  Ag°T  j 
i 


(10) 
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where  R  is  the  universal  gas  constant,  AG°  is  the  standard  Gibbs 
function  of  reaction  and  A g°f  Ti  represents  the  standard  Gibbs 
function  of  formation  at  given  temperature  T  (K)  of  the  gas  species 
i  which  can  be  expressed  by  the  empirical  equation  given  below: 

A  g°fj  =  hf  -  aT  In  (T)  -  bT2 -^T3 -^T4+^f+f+gT  (11) 

The  values  of  coefficients  a-g  and  the  enthalpy  of  formation 
(hj)  of  different  gases  are  available  in  the  literature  [9], 


initial  elemental  abundance  vector  and  its  current  iteration  value 
(bk-b[m\  is  added  to  the  right  hand  side  of  Eq.  (15)  to  eliminate 
error  accumulation  during  the  iteration  process  [14]. 

Finally,  the  new  numbers  of  moles  vector,  n(m+1),  is  determined  by 

n(m+l)  =  n(m)+(B(m)5n(m)  (18) 

where  ®(m)  is  the  step  size  parameter. 

Finally,  the  molar  fractions  of  all  species  are  determined  by 


3.1.2.  Non-stoichiometric  equilibrium  model 

The  non-stoichiometric  equilibrium  model  is  based  on  the 
minimization  of  Gibbs  free  energy  of  the  system  without  specifi¬ 
cation  of  the  possible  reactions  taking  place.  Here,  only  the 
elemental  composition  of  the  fuel  is  needed  to  specify  which  can 
be  obtained  from  the  ultimate  analysis  of  the  fuel.  This  method  is 
particularly  suitable  for  problems  with  unclear  reaction  mechan¬ 
isms.  A  non-stoichiometric  equilibrium  modeling  approach  based 
on  the  literature  [10,11]  is  presented  below.  The  model  is  based  on 
the  RAND  [12,13]  algorithm.  Here,  the  change  in  mole  of  a  species 
in  the  mth  iteration  can  be  expressed  explicitly  as  a  function  of  its 
current  chemical  potential,  the  phase  distribution  of  the  species  at 
a  given  system  temperature  and  pressure  and  the  Fagrange 
multiplier  as  given  below: 

for  multi  -  species  phase 

=  u„njm)  for  single  -  species  phase 

(1  =  1, 2,  ...,N;  fc=  1,2,  ...,/C;  a=  1, 2,  ...,*)  (12) 

Here,  N  denotes  the  total  number  of  species,  n(m)  denotes  the  mole 
of  species  i  in  the  inth  iteration,  aik  is  the  coefficient  in  species- 
element  matrix,  ^-m)  is  the  chemical  potential  of  species  i  in  the 
mth  iteration,  i//k  is  a  function  related  to  Lagrange  multiplier  Ak  and 
ua  is  the  phase  split  of  <5njm). 

Here,  i//k  is  given  by 

Vk  =  ^f  (13) 

and  ua  is  given  by 

ua=  2  sn™/nr  =  dn(£)/nr  (14) 

i  =  1 


(5n<m»  =  n!m) 


* 

2  aik¥k~l~^a  nT 
k  =  1 


where  subscript  f  means  total  and  a  refers  to  the  phase  to  which  a 
species  belongs. 

Therefore,  a  set  of  (K+x)  simultaneous  algebraic  equations  that 
are  to  be  solved  iteratively  by  the  RAND  algorithm  includes  K 
linear  equations  regarding  element  abundance 


2  2  a*a„n?"yk+  f  bj>, 

k=\i=\  a  =  1 

0  =  1,2,. ..,/C) 


N 

Z  aijJ 

i  =  1 


F%r+bk-tiP 


RT 


(15) 


and  7i  supplementary  equations  for  different  phases 


2  b™Vk-nzau, 

k=  1 

(a=  1,2,  ...,x) 


N  „<m> 

y  n(m)f‘ia 

ifi  RT 


(16) 


The  initial  elemental  abundance  vector  b  is  calculated  from  the 
feedstock  and  the  fcth  element  of  the  b  vector  at  the  mth  iteration  is 

C=  2  (17) 

i  =  1 

Mass  balance  constraints  are  imposed  at  every  iterations 
during  solution  of  above  equations,  while  the  algorithm  iteratively 
minimizes  the  Gibbs  free  energy.  The  difference  between  the 


Equilibrium  models  are  extensively  used  by  many  authors  [15-23] 
due  to  its  advantageous  feature  of  simplicity  and  it  gives  quick  idea 
of  limits  of  operation  compared  to  the  rate  base  model.  It  can 
provide  the  final  gas  composition  variation  with  operating  condi¬ 
tion,  gas  yield,  calorific  value  of  gas,  carbon  conversion  etc.  The 
prediction  using  equilibrium  model  matches  well  with  the  experi¬ 
ment  when  the  equilibrium  condition  is  prevailed  i.e.  when  the  rate 
of  reaction  is  very  fast  compared  to  the  residence  time  of  reactants 
in  the  gasifier.  The  equilibrium  model  cannot  predict  the  spatial 
variation  of  the  product  inside  the  gasifier  as  well  as  the  temporal 
variation.  It  is  also  important  to  note  that  the  predicted  methane 
concentration  from  equilibrium  model  is  always  substantially  lower 
than  the  experimentally  observed  data  for  gasification  process. 
The  high  amount  of  methane  in  the  experimental  gas  composition 
is  a  non-equilibrium  species  resulting  from  pyrolysis  reactions  and 
limited  char  gasification  is  achieved  in  the  gasifier.  So,  the  measured 
gas  composition  is  largely  attributed  to  pyrolysis  which  the  equili¬ 
brium  model  cannot  account  for.  The  equilibrium  model  cannot 
predict  highly  accurate  result  for  fluidized  bed  gasifier  due  to  its 
inherent  assumptions.  The  residence  time,  feed  configuration, 
internal  mixing  and  gas-solid  contacting  pattern  in  fluidized  bed 
gasifier  account  for  large  deviation. 

Thus,  in  order  to  predict  more  accurately  the  behavior  of 
fluidized  bed  gasifier,  different  modifications  of  the  existing 
equilibrium  models  have  been  done  as  discussed  below.  Li  et  al. 
[10,11]  modified  their  non-stoichiometric  equilibrium  model  by 
introducing  the  experimental  carbon  conversion  from  fluidized 
bed  gasification  and  got  better  agreement  to  the  experimental 
data.  For  example,  they  compared  the  higher  heating  value  (HHV) 
of  the  product  gas  as  a  function  of  air  ratio  from  experiment,  pure 
equilibrium  model  and  modified  equilibrium  model  as  shown 
in  Fig.  3.  It  was  observed  that  the  modified  equilibrium  model 


Fig.  3.  Effect  of  air  ratio  on  dry  gas  HHV  for  as-received  high  value  coal  gasified  at 
155  kPa  and  1020-1150  K.  Points  are  experimental  value  [9]. 
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Table  1 

Gas  composition  by  using  pure  equilibrium  model  and  modified  equilibrium  model  and  experiment  [5], 


Tempe  (°C)  S/B  Experiment  Pure  equilibrium  model  RMS  error  Modified  equilibrium  model  RMS  error 


h2 

CO 

co2 

ch4 

h2 

CO 

n 

o 

ch4 

h2 

CO 

co2 

ch4 

690 

1.32 

50.50 

14.30 

26.60 

8.60 

48.76 

12.99 

31.21 

7.03 

2.667429 

48.4 

16.26 

28.88 

6.46 

2.123064 

730 

1.32 

52.20 

16.40 

23.50 

7.90 

50.71 

15.48 

28.74 

5.06 

3.10603 

50.14 

18.94 

26.33 

4.6 

2.720023 

750 

1.00 

49.50 

23.70 

21.20 

5.60 

50.37 

20.59 

25.01 

4.02 

2.619232 

49.64 

24.22 

22.52 

3.62 

1.219918 

750 

1.32 

52.30 

17.75 

22.25 

7.40 

51.43 

16.64 

27.65 

4.28 

3.197006 

50.76 

20.15 

25.22 

3.87 

2.711706 

750 

1.70 

52.90 

16.40 

22.90 

7.80 

52.26 

13.50 

29.74 

4.50 

4.077229 

51.67 

16.78 

27.47 

4.08 

3.015817 

770 

1.32 

54.40 

18.50 

19.40 

7.70 

52.00 

17.72 

26.66 

3.62 

4.350931 

51.24 

21.28 

24.23 

3.25 

3.900173 

Average  RMS  3.33631  Average  RMS  2.615117 


predicted  the  experimental  values  much  better  than  the  pure 
equilibrium  model.  Loha  et  al.  [6]  studied  the  fluidized  bed 
biomass  gasification  by  using  a  stoichiometric  equilibrium  model. 
The  pure  equilibrium  model  was  modified  by  multiplying  equili¬ 
brium  constants  with  pre  factors.  Results  from  the  pure  equili¬ 
brium  model,  the  modified  equilibrium  model  and  the  experiment 
were  compared  as  shown  in  Table  1.  It  was  observed  that  average 
RMS  error  reduced  from  3.33631  to  2.615117  while  using  the 
modified  equilibrium  model.  Fryda  et  al.  [24]  and  Schuster  et  al.  [25] 
also  evaluated  the  product  gas  composition  from  fluidized  bed 
gasification  by  using  equilibrium  model.  They  modified  the  equili¬ 
brium  model  by  introducing  the  un-reacted  char  carbon  as  some 
percentage  of  the  biomass  carbon  input  and  got  better  agreement  to 
the  experimental  data. 

Therefore,  if  the  equilibrium  model  is  modified,  it  can  be  use  to 
study  the  fluidized  bed  gasification  performance,  particularly  to 
study  the  influence  of  most  important  parameters  on  the  overall 
performance.  Equilibrium  model  also  provides  the  thermodynamic 
limits  of  the  gasification  system.  Hence,  the  equilibrium  model 
may  be  useful  for  preliminary  design  and  comparison  of  fluidized 
bed  gasification  process  if  they  are  modified  as  mentioned  above. 
But,  it  cannot  provide  highly  accurate  data  for  all  cases  because  it 
does  not  consider  the  hydrodynamic  behavior  of  the  system  as 
well  as  the  reaction  kinetics.  Hence,  to  get  more  accurate  data  for 
detail  design,  the  rate  base  modeling  or  kinetic  modeling  with 
hydrodynamic  consideration  of  the  fluidized  bed  system  is 
required  as  discussed  in  the  following  sections. 

3.2.  Rate  base  model  or  kinetic  model 

Rate  base  model  or  kinetic  model  is  more  comprehensive  and 
realistic  compared  to  the  equilibrium  model.  It  takes  into  account 
the  hydrodynamics,  transport  process  and  reaction  kinetics.  In  rate 
base  model,  a  detailed  reaction  mechanism  with  rate  of  each 
reactions  are  considered  depending  on  the  type  of  fuel  used. 
Modeling  the  reaction  kinetics  in  gasification  process  requires 
the  consideration  of  complex  reaction  network  which  consists  of 
pyrolysis,  heterogeneous  gas-solid  reactions  and  homogeneous 
gas  phase  reactions  as  discussed  below. 

3.2.1.  Pyrolysis 

The  solid  fuel  particles  undergo  pyrolysis  upon  entering  into 
the  hot  fluidized  bed  gasifier.  Pyrolysis  is  basically  the  decomposi¬ 
tion  of  solid  fuel  resulting  from  heat.  In  gasification,  the  final 
product  gas  composition  is  important  because  it  determines  the 
heating  value  of  the  product  gas.  Therefore,  the  kinetic  model  for 
total  devolatilization,  which  determines  both  the  extent  of  pro¬ 
duced  char  and  individual  char  released,  should  be  taken  into 
consideration. 

Kinetic  models  used  for  determining  the  pyrolysis  product, 
available  in  the  literature,  are  either  based  on  a  single  equation 


model,  two-equation  model  or  a  combination  of  reactions.  The  single 
equation  model  is  the  simplest  one  and  it  is  represented  as  [26] 


Fueli//  volatile +(1  ->i)  char  (20) 

k=Aexp(m)  (21) 

In  two  equations  model,  the  pyrolysis  of  any  solid  fuel  is  represented 
by  two  equations  as  given  below  [27]: 

Fuel^i;/!  volatile + (1 -?/])  char  (22) 

Fuel^;?2  volatile + (1  — >/2)  char  (23) 

where 

ki=AeXp  (etJ)  (24) 


Parameters  Land  >/  in  the  above  equations  are  to  be  obtained 
from  experiments.  Many  publications  deal  with  the  calculation  of 
those  parameters  for  various  fuels,  but  there  is  a  great  variation 
between  the  representations,  probably  because  a  number  of 
physical  and  chemical  factors  are  incorporated  in  one  expression. 
Despite  this  fact,  this  type  of  expression  has  been  widely  used  in 
reactor  models  due  to  its  simplicity  making  it  computationally 
tractable,  needing  only  a  small  set  of  input  data.  In  most  cases, 
however,  it  does  not  give  enough  information  for  comprehensive 
modeling. 

The  models  mentioned  above  can  predict  the  amounts  of  volatile 
and  char  released  during  pyrolysis,  but  the  yields  of  the  main  gas 
species  are  not  predicted  which  also  varies  with  the  type  of  fuel. 
The  pyrolysis  gas  composition  is  generally  modeled  either  based  on 
the  experimental  data  available  from  the  literature  for  a  particular 
fuel  [28-30]  or  based  on  the  developed  correlation  [31-33],  These 
correlations  are  also  developed  for  a  particular  set  of  operating 
conditions  and  type  of  fuel  which  is  not  always  similar  to  those  has 
to  be  simulated.  Therefore,  more  fundamental  data  are  desirable  to 
increase  the  predictive  capability  of  model  simulations. 

3.2.2.  Gasification  with  air 

In  case  of  gasification  with  air,  the  heterogeneous  reactions  of 
carbon  with  oxygen  in  the  air  are  highly  exothermic  and  therefore 
provide  the  heat  required  for  the  subsequent  gasification  reac¬ 
tions.  The  heterogeneous  reactions  of  carbon  with  oxygen  are 
given  below 

C+02^C02,  Ah°  =  -393.5  kj/mol  (25) 

C+^02^C0,  Ah°R=  -110.5  kj/mol  (26) 

In  addition  to  the  heterogeneous  gas-solid  reactions,  there  are  also 
many  homogeneous  gas  phase  reactions.  These  are  the  combustion 
reactions  of  all  the  gases  produced  during  the  pyrolysis  step. 
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The  hydrogen  and  carbon  monoxide  will  be  oxidized  completely  to 
carbon  dioxide  and  water  as  given  below: 

H2+^02^H20,  Ah°R  =  -241.8  kj/mol  (27) 

CO+lo2^C02,  AhR  =  -283.0  kj/mol  (28) 

In  addition,  the  complete  combustion  of  hydrocarbons  with  oxygen 
will  produce  carbon  dioxide  and  water.  But,  in  the  reducing  atmo¬ 
sphere  of  gasifier  and  lack  of  oxygen,  the  producs  will  be  hydrogen 
and  carbon  monoxide. 

CH4  +  ^02-C0  +  2H2,  AhR  =  -35.7  kj/mol  (29) 

C2H4  +  02  ->2CO+2H2,  Ah°  =  -35.7  kj/mol  (30) 

In  gasification,  the  amount  of  air  added  to  the  gasifier  is  not 
sufficient  to  completely  combust  the  fuel  to  carbon  dioxide  and 
water.  Therefore,  once  all  oxygen  is  consumed,  further  reaction 
taking  place  are  homogeneous  equilibrium  reactions  or  heteroge¬ 
neous  reactions  between  char  carbon  and  produced  gases  so  far. 
Therefore,  the  gasification  with  air  is  basically  a  combination  of 
steam  and  carbon  dioxide  gasification. 

3.2.3.  Gasification  with  steam 

If  steam  is  used  as  the  gasifying  agent  or  all  oxygen  in  air 
gasification  is  consumed,  the  gasification  reactions  with  steam  are 
the  main  reactions  in  the  temperature  range  normally  used  in 
fluidized  bed.  The  products  of  the  reactions  of  solid  carbon  with 
steam  are  hydrogen,  carbon  monoxide  and  to  a  smaller  extent  also 


carbon  dioxide  may  be  produced. 

C+H2O^CO+H2,  Aft°=  +131.3  kj/mol  (31) 

C  +  2H20->C02  +  2H2,  Ah°R  =  +211.4  kj/mol  (32) 

The  homogeneous  gas  phase  reaction  of  steam  with  methane  is 
given  below: 

CH4  +  H2O^CO  +  3H2  AhR  =  +206.2  kj/mol  (33) 

The  very  well  known  water-gas-shift  reaction  which  is  indepen¬ 
dent  of  the  choice  of  the  gasifying  agent  and  operating  condition  is 
given  below: 

C0+H20<-»C02  +  H2,  AhR  =  -41.2  kj/mol  (34) 


Other  homogeneous  steam  gasification  reactions  are  the  reactions 
of  steam  with  higher  hydrocarbons.  The  main  product  of  these 
reactions  will  be  other  hydrocarbons,  carbon  monoxide  and 
hydrogen. 

3.2.4.  Gasification  with  carbon  dioxide 

The  gasification  reaction  with  carbon  dioxide  requires  higher 
temperature  because  this  reaction  proceeds  veiy  slowly.  The 
heterogeneous  gasification  reaction  of  carbon  dioxide  with  char 
carbon  is  given  below: 

C+C02 -»2CO,  Ah°R=  +172.5  kj/mol  (35) 

The  homogeneous  gas  phase  reaction  of  methane  with  carbon 
dioxide  is  given  below. 

CH4  +  C02^2C0  +  2H2,  AhR=  +247.3  kj/mol  (36) 

3.2.5.  Other  gasification  reactions 

In  addition  to  the  reactions  discussed  above,  there  may  be 
some  other  reactions  in  the  gasification  process.  For  example,  the 
heterogeneous  gasification  reaction  of  char  carbon  with  hydrogen 
is  given  below.  This  reaction  is  significant  only  at  higher  pressure: 

C+2H2->CH4,  AhR  =  —74.9  kj/mol  (37) 


4.  Fluid  dynamics 

Hydrodynamics  plays  an  important  role  in  defining  the  perfor¬ 
mance  of  fluidized  bed  gasifier.  Fluidized  bed  exhibits  very  com¬ 
plex  hydrodynamics  due  to  the  nonlinear  interactions  between 
fluid  and  particle  and  their  own  individual  movement  tendencies. 
Depending  on  the  description  of  the  hydrodynamics,  the  fluidized 
bed  gasification  models  may  be  classified  as  the  two-phase  flow 
model,  the  Euler-Euler  model  and  the  Euler-Lagrange  model. 
Detailed  descriptions  of  each  of  these  models  are  presented  below. 

4.3.  Two-phase  flow  model 

Two  most  popular  and  oldest  two-phase  flow  models  are  as  pro¬ 
posed  by  Davidson  and  Harrison  [34]  and  Kunii  and  Levenspiel  [35], 
In  two-phase  flow  model,  the  fluidized  bed  is  divided  into  two 
phase  e.g.  the  bubble  phase  and  the  emulsion  phase  (Fig.  4).  The 
bubble  phase  is  assumed  to  be  particle  free.  The  emulsion  phase 


Fig.  4.  Division  of  fluidized  bed  into  bubble  and  emulsion  phases  in  two-phase  flow  model. 
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contains  particles  and  remains  in  minimum  fluidization  condition. 
Gas  flows  in  plug  flow  in  both  the  bubble  and  emulsion  phase.  The 
solids  are  evenly  distributed  and  are  transported  by  dispersion  only. 
The  heterogeneous  gas-solid  reactions  occur  in  the  emulsion  phase 
only.  Homogeneous  gas  reactions  take  place  in  both  the  bubble  and 
the  emulsion  phase.  The  model  posses  separate  mass  balance 
equation  for  each  phase  and  mass  interaction  between  phases.  In 
two-phase  flow  model,  the  momentum  equation  is  not  explicitly 
solved  rather  the  hydrodynamic  of  the  fluidized  bed  is  calculated 
from  correlations.  While  modeling  fluidized  bed  by  using  two-fluid 
model,  the  gasifier  is  axially  divided  into  different  zones  (bottom 
zone,  freeboard  and  exit  zone).  The  freeboard  and  the  exit  zone  are 
generally  modeled  by  comparatively  simple  mixture  model.  The  two- 
phase  flow  concept  is  used  to  model  the  dense  bottom  zone  in 
fluidized  bed  systems.  Mathematical  equations  used  for  such  a  two- 
phase  flow  model  following  the  work  of  Kruse  and  Werther  [36]  and 
Petersen  and  Werther  [22,37]  are  given  below. 

In  two-phase  flow  model,  separate  mass  conservation  equa¬ 
tions  are  solved  for  gas  and  solid  species.  Gaseous  species  are 
present  both  in  bubble  and  emulsion  phases  and  hence  separate 
mass  conservation  equations  are  solved  for  gas  species  in  bubble 
and  emulsion  phases. 

The  conservation  equation  for  the  gaseous  component  fin  the 
bubble  phase  is  given  by 

a  t)C^  +  ()(uhChj)  =  - l<bACghj _ Cged + ^  2  v..rbj  (38> 

01  OL  jig-g ) 

The  conservation  equation  for  the  gaseous  component  i  in  the 
emulsion  phase  is  given  by 

,  dCze  i  .  d(Cse  j)  d2(Cse ,) 

(1  +  -ab)Umf— — - (1  -ab)amfDgtax  ^  ' 

=  Kbe(Cgb,i  —  Cge  i)  -{-  (1  —  ab)ocmf  2  ^ijl~emulj~^~6v  bz  ^  V ijTemulJ 

j(g-s)  iig-g ) 

(39) 


The  particles  are  present  in  the  emulsion  phase  only.  Hence,  the 
mass  balance  equation  of  solid  component  /in  the  emulsion  phase 
is  given  by 


oxs,i  n  °  (*s,i 

LvPs  LvPsUS,OX  ^2 


where  ab  is  the  void  fraction  of  bubble  phase,  amt  is  the  void 
fraction  at  minimum  fluidization  condition,  cv  is  the  solid  volume 
concentration,  ub  is  the  velocity  of  gas  in  the  bubble  phase,  um!  is 
the  minimum  fluidization  velocity,  Cgb  i  is  the  concentration  of  gas 
species  /  in  the  bubble  phase,  Cgei  is  the  concentration  of  gas 
species  i  in  the  emulsion  phase,  xsi  is  the  mass  fraction  of  solid 
species  i,  Kbe  is  the  mass  transfer  coefficient  between  bubble  and 
emulsion  phase,  vy  is  the  stoichiometric  coefficient  of  species  i  in 
reaction  j,  Dgax  is  the  axial  dispersion  coefficient  of  gas  phase, 
Dsax  is  the  axial  dispersion  coefficient  of  solid,  rbj  is  the  reaction 


rate  of  reaction  jin  bubble  phase,  remuy  is  the  reaction  rate  of 
reaction  j  in  emulsion  phase,  rh'"Feedi  is  the  volumetric  mass  flow 
rate  at  the  feeding  section,  m'"Retumi  is  the  volumetric  mass  flow 
rate  at  the  return  leg  and  M,  is  the  molar  mass  of  species  /. 

The  above  mass  balance  equations  are  discrelized  and  numeri¬ 
cally  solved  with  appropriate  initial  and  boundary  conditions. 
In  two-phase  flow  model,  the  fluid  dynamic  parameters  are 
calculated  based  on  correlations  and  supplied  as  input  parameters. 
Several  expressions  exist  to  calculate  these  parameters.  A  list  of 
hydrodynamic  parameters  generally  used  in  two-phase  flow 
modeling  is  given  in  Table  2.  Therefore,  the  need  to  solve  the 
momentum  equations  is  avoided  and  the  model  becomes  com¬ 
paratively  simple.  A  typical  solution  procedure  for  two-phase  flow 
modeling  of  fluidized  bed  gasification  process  is  shown  in  Fig.  5. 

Several  studies  are  reported  on  fluidized  bed  gasification  by 
using  the  two-phase  flow  modeling  concept  discussed  above  or  in 
modified  form  by  assuming  or  estimating  several  terms.  Jiang  and 
Morey  [33]  studied  the  fluidized  bed  biomass  gasification  by  using 
one  dimensional  steady-state  two-phase  flow.  The  model  pre¬ 
dicted  fuel  feed  rate,  composition  of  product  gas  and  fuel  energy 
conversion.  Model  results  were  compared  with  experimental  date. 
It  was  observed  that  model  worked  well  at  higher  temperature 
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i 
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Fig.  5.  Flowchart  for  solution  procedure  of  two-phase  flow  model. 


Table  2 

Correlations  used  for  calculating  the  hydrodynamic  parameters  in  two-phase  flow  model. 


A  widely  used  correlation  given  by  Wen  and  Yu  [38]  for  calculating  umfis  given  by 

umf  =  Re’’  mf/Jg ,  Where  Repm/  =  <Jc\+C2At - C,  and  Ar  =  i’g>S 

aPPg 

Several  values  for  C\ and  C2  are  proposed  in  the  literature  [39]. 

The  void  fraction  at  minimum  fluidization  ( amf )  can  be  calculated  by  solving  the  equation  derived  from  Ergun  equation  [79] 
!-75  p  2  150(1  -amf) 

w  Re”"lf-Ar 

The  void  fraction  of  bubble  ( ab )  can  be  calculated  as 


<*b  =  1  - 


( 1  amf ) 


The  solid  volume  concentration  cv  can  be  calculated  as 

Cv  =  (1  —  «b)(l  —  amf) 
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Fig.  6.  Effect  of  oxygen  in  the  feed  on  (A)  bed  temperature  and  calorific  value  and 
(B)  product  gas  composition.  Steam:  6.8  kg/h;  Coal:  coke-breeze,  10  kg/h  [40], 

and  for  low  temperature  it  was  not  good.  Chatterjee  et  al.  [40] 
studied  the  gasification  of  high  ash  Indian  coal  in  a  lab-scale 
fluidized  bed  gasifier  by  using  two-phase  flow  model.  Steam-and 
air  was  used  as  the  gasifying  agent.  The  gas  composition,  tem¬ 
perature,  carbon  conversion  and  calorific  value  of  gas  with 
different  operating  conditions  were  simulated.  The  comparison 
between  simulated  and  experimental  results  showed  similar 
trend.  For  example,  the  simulated  and  experimental  gas  composi¬ 
tion,  calorific  value  and  temperature  as  a  function  of  oxygen  feed 
rate  for  coke-breeze  and  bituminous  coal  were  shown  (Fig.  6)  and 
the  estimated  error  between  prediction  and  experiment  was 
within  12%.  Yan  et  al.  [41]  studied  the  fluidized  bed  coal  gasifica¬ 
tion  by  using  two-phase  flow  model  and  showed  that  the  change 
in  volume  of  gas  due  to  reaction  resulted  in  a  higher  velocity.  The 
suspension  phase  gas  stayed  at  minimum  fluidization,  only  the 
bubble  phase  gas  velocity  was  increased.  Thus,  a  net-flow  concept 
was  accounted  for  which  considered  the  production  or  reduction 
of  gaseous  volume  by  reaction  and  it  was  directed  from  the 
suspension  phase  to  the  bubble  phase.  The  results  showed  that 
the  net-flow  is  significant,  in  the  range  71-87  %,  relative  to  the 
feed  gas  rate.  Simulation  without  net-flow  deviated  significantly 
from  experiment.  Jennen  et  al.  [42]  studied  the  gasification  of 
wood  in  a  pilot-scale  circulating  fluidized  bed  by  using  one¬ 
dimensional  two-phase  flow  model.  The  calculated  axial  profile 
of  the  gas  composition  and  the  temperature  are  compared  with 
the  experimental  data  and  very  good  match  was  found.  The 
difference  between  the  measured  and  the  calculated  gas  composition 


was  +1%  and  maximum  deviation  of  temperature  was  5  °C.  Hamel 
and  Krumm  [43]  used  the  two-phase  flow  model  to  simulate  four 
bubbling  fluidized  bed  gasifiers  of  different  scale  from  atmospheric 
lab-scale  to  pressurized  commercial-scale.  Simulated  results  of  over¬ 
all  carbon  conversion,  freeboard  temperature  and  the  gas  species 
concentration  deviated  + 10%  from  experiment.  Fiaschi  and  Miche- 
lini  [44]  used  a  two-phase  flow  model  to  study  the  biomass 
gasification  kinetics  in  bubbling  fluidized  bed  gasifier.  The  model 
predicted  temperature  and  concentration  gradient  along  the  axis. 
The  model  result  showed  largely  satisfactory  agreement  with  the 
experimental  result  from  the  literature  and  the  requirement  of 
further  validation  of  the  model  was  identified.  Sadaka  et  al.  [45-47] 
developed  a  two-phase  model  for  predicting  the  performance  of 
dual  distributor  type  fluidized  bed  biomass  gasifier.  They  assumed 
that  the  fluidized  bed  consisted  of  a  dilute  phase  (jets,  bubbles 
and/or  slugs)  and  an  emulsion  phase.  The  emulsion  phase  was 
divided  into  an  interstitial  gas  phase  and  a  solid  phase.  Model 
predicted  the  bed  temperature,  gas  composition,  heating  value 
and  gas  yield.  Sensitivity  analysis  was  carried  out  by  varying 
fluidization  velocity,  steam  flow  rate  and  biomass-to-steam  ratio. 
The  model  results  predicted  the  experimental  data  with  higher 
accuracy  (R2= 0.88  -0.98).  Chejne  and  Hernandez  [48]  developed 
a  one  dimensional  two-phase  model  to  simulate  the  coal  gasifica¬ 
tion  in  fluidized  bed.  The  model  predicted  temperature,  gas 
composition,  volume  fraction,  velocity  and  other  fluid  dynamic 
parameters.  The  model  results  showed  20%  error  while  compared 
with  the  experimental  results  from  the  literature.  Ross  et  al.  [49] 
improved  the  isothermal  two-phase  flow  model  by  considering 
the  non-isothermal  behavior  of  gases  and  heat  transfer  mechan¬ 
isms  and  used  to  study  the  fluidized-bed  coal  gasification  process. 
Isothermal  and  non-isothermal  simulations  were  run  for  the  pilot- 
scale  and  commercial-scale  gasifier.  Gas  composition,  bed  tem¬ 
perature  and  reaction  rate  were  predicted.  The  good  agreement 
with  experimental  results  was  found  for  non-isothermal  simula¬ 
tion  and  significant  deviation  was  found  for  isothermal  simulation. 
Petersen  and  Werther  [22]  used  two-phase  flow  models  for 
simulating  the  gasification  of  sewage  sludge  in  a  pilot-scale 
circulating  fluidized  bed  gasifier.  Initially,  a  pseudo-two- 
dimensional  model  was  developed  and  the  gas  composition  was 
calculated  by  selecting  the  reaction  kinetics  from  the  literature. 
Then  the  kinetics  of  pyrolysis  and  main  gasification  reactions  were 
determined  by  comparing  the  pseudo-two-dimensional  model 
results  with  the  experimental  data.  Finally,  a  three-dimensional 
model  [36]  was  developed  with  the  modified  reaction  kinetics. 
Fig.  7  shows  the  axial  gas  composition  profile  for  pyrolysis,  C02 
gasification  and  air-gasification  obtained  from  experiment  and 
pseudo-two-dimensional  model  with  fitting  parameters.  Good 
agreement  between  calculated  and  measured  gas  composition 
was  seen  after  fine  tuning  the  reaction  kinetics.  Radmanesh 
et  al.  [50]  used  one  dimensional  isothermal  two-phase  flow  model 
to  study  the  fluidized  bed  gasification  of  biomass.  Two  different 
kinetic  models  for  pyrolysis  were  used  and  their  impact  was 
studied.  The  effect  of  equivalence  ratio,  steam-to-biomass  ratio, 
bed  temperature,  feed  location,  and  mass  transfer  between  the 
countercurrent  regions  on  the  gas  composition  was  studied.  They 
showed  the  importance  of  the  pyrolysis  step  in  predicting  the  final 
gas  composition  by  comparing  the  simulated  gas  composition  with 
experimental  data.  The  pyrolysis  model  derived  at  higher  heating 
rates  estimated  the  final  gas  composition  relatively  better  than  the 
one  derived  at  lower  heating  rate.  The  strong  influence  of  pyrolysis 
step  on  the  overall  performance  of  the  fluidized  bed  biomass 
gasification  was  also  identified  by  Kaushal  et  al.  [51]  in  their  two- 
phase  flow  modeling  study.  Pengmei  et  al.  [52]  studied  the 
biomass  gasification  in  fluidized  bed  using  two-phase  flow  model 
and  showed  that  the  trend  of  changing  the  gas  composition  with 
temperature  was  in  accordance  with  the  experiment.  However,  the 
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Fig.  7.  Axial  profile  of  simulated  and  measured  gas  composition  in  (A)  pyrolysis,  (B)  C02  gasification,  and  (C)  air-gasification  experiments  and  modeling  results  from  pseudo- 
two-dimensional  model  [36], 


gas  yield  from  the  simulation  was  higher  that  experiment.  It  was 
explained  that,  in  simulation,  all  char  particles  were  assumed  to 
take  part  in  the  reaction  and  no  illusion  or  escape  from  the  bed 
which  resulted  more  gas  yield.  Goyal  et  al.  [53]  studied  the 
fluidized  bed  gasification  of  mixture  of  coal  and  petcoke  by  using 
two-phase  flow  model.  The  effect  of  various  operating  parameters 
such  as  composition  of  feed,  location  of  feed  point  and  ash  content 
on  the  performance  of  the  gasifier  were  studied.  Results  showed 
that  the  increase  in  petcoke  content  in  the  feed  mixture  tends  to 
lower  the  efficiency  and  carbon  conversion  but  increases  the 
amount  of  syngas  produced  and  increase  in  ash  content  of  coal 
decreases  the  carbon  conversion.  The  model  results  also  identified 
that  feed  point  of  the  solids  should  be  above  the  point  where  02  in 
the  bed  gets  exhausted,  in  order  to  obtain  the  maximum  carbon 
conversion  and  efficiency.  Gungor  [54]  used  two-phase  flow 


model  to  study  the  fluidized  bed  biomass  gasification  process 
including  the  tar  conversion.  Effects  of  gasification  temperature,  bed 
operational  velocity,  equivalence  ratio,  biomass  particle  size  and 
biomass-to-steam  ratio  on  the  hydrogen  production  was  studied. 
The  model  was  capable  of  predicting  concentration  distribution  of  gas 
species  and  tar  along  the  height  of  the  gasifier.  Simulated  hydrogen 
production  was  compared  with  experimental  data  which  showed 
largely  satisfactory  agreement  but  the  effect  of  considering  the  tar 
conversion  was  not  shown. 

Till  date,  the  two-phase  model  is  the  most  established  model 
for  simulating  the  fluidized  bed  gasification  process  published  in 
the  open  literature.  It  is  comparatively  simple  because  it  does  not 
consider  the  complex  gas-solid  dynamics  but  still  maintains  the 
fluid  dynamic  effect.  Therefore,  it  can  provide  important  informa¬ 
tion  about  the  fluidized  bed  gasifier  like  the  transient  nature  of  the 
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gasifier,  distribution  of  gas  species  and  reacting  particles  inside  the 
reactor,  bubble  size  and  velocity,  temperature  distribution  etc.  But, 
due  to  the  fact  that  the  hydrodynamics  is  based  on  some  correla¬ 
tions,  the  applicability  of  the  two-fluid  model  is  restricted  to  the 
operating  condition  on  which  these  correlations  are  valid.  The 
most  common  assumption  of  two-phase  flow  model  is  the  perfect 
mixing  of  solid  and  plug  flow  of  gas  in  bubbles  and  emulsion 
phases  which  are  sometimes  not  prevail  in  the  bed.  The  two-phase 
flow  modeling  approach  cannot  be  applied  to  the  freeboard  region 
of  the  bed,  the  freeboard  is  separately  modeled.  Besides,  most  of 
the  two-phase  flow  models  are  either  one  dimensional.  Therefore, 
the  computational  fluid  dynamic  (CFD)  model  by  coupling  the 
detail  fluidized  bed  hydrodynamics  and  reaction  kinetics  may 
overcome  the  above  disadvantages  of  two-fluid  model  and  can 
provides  more  inside  into  the  fluidized  bed  gasifier.  In  CFD  models, 
the  hydrodynamics  of  fluidized  bed  is  based  on  either  the  Euler- 
Euler  concept  or  the  Euler-Lagrange  concept  as  discussed  below. 

4.2.  Euler-Euler  model 

In  Euler-Euler  model,  the  gas  and  the  solid  phases  are  treated 
as  interpenetrating  continuum  present  at  the  same  time  in  the 
same  control  volume.  The  hydrodynamic  model  uses  the  conser¬ 
vation  of  mass  and  momentum  of  each  phase  and  represented  by 
respective  conservation  equations.  The  interaction  between  two 
phases  is  expressed  as  additional  source  terms  added  to  the 
conservation  equations.  The  conservation  equations  for  mass  and 
momentum  are  solved  for  the  gas  and  solid  phase  separately  and 
are  closed  with  appropriate  constitutive  relations. 

Continuity  equation  for  the  gas  phase  is  given  by 

— (cigPg)  +  V.(agpg  v  g)  =  SgS  (41 ) 

Continuity  equation  for  the  solid  phase  is  given  by 

^(“sPs)  +  V.(asps"v  s)  =  Sgs  (42) 

where  a  is  the  volume  fraction,  p  is  the  density,  v  is  the  velocity 
vector  and  Sgs  is  the  mass  source  term  due  to  chemical  reactions. 
The  momentum  equation  for  gas  phase  is  given  by  the  Navier- 
Stokes  equation,  modified  to  include  the  inter-phase  momentum 
transfer  term. 

Momentum  equation  for  the  gas  phase  is  given  by 

~^(agPg  v  g )  +  V. («„/)„  v  g  v  g)  =  —  agVp  +  V.Tg+agpg  g  +Kgs(  vs-  v  g)-t-Sgs  v  g 

(43) 

where  p  is  the  hydrodynamics  pressure,  zg  is  the  viscous  stress 
tensor,  g  is  the  acceleration  due  to  gravity  and  Kgs  is  the 
momentum  transfer  between  gas  and  solid  phases. 

The  stress  tensor  of  the  gas  phase  can  be  given  as 

.  2 

Tg  =  agPg(^  v  g  +  ^  v  g)  —  j  “g  Eg(y  ■  v  g)^  (44) 


where  ps  is  the  solid  shear  viscosity  and  Xs  is  the  solid  bulk 
viscosity. 

The  Euler-Euler  model  became  more  popular  after  the  devel¬ 
opment  of  Kinetic  Theory  of  Granular  Flow  (KTGF)  which  is  based 
on  the  theory  of  non-uniform  dense  gases  described  in  the 
literature  [55].  The  pioneering  paper  of  Lun  et  al.  [56]  applied 
the  kinetic  theory  of  gases  to  granular  flow.  Sinclair  and  Jackson  [57] 
applied  the  granular  flow  model  to  a  fully  developed  gas-solids  flow 
in  a  pipe.  Ding  and  Gidaspow  [58]  derived  expressions  for  solids 
viscosity  and  pressure  of  a  dense  gas-solids  flow.  Gidaspow  [59] 
extended  the  Ding  and  Gidaspow  [58]  formulation  to  both  dilute 
and  dense  phase  by  considering  a  non-Maxwellian  velocity  distribu¬ 
tion.  Huilin  et  al.  [60]  applied  KTGF  to  study  the  motion  of  particles 
in  the  gas  bubbling  fluidized  bed  with  the  binary  mixtures.  Louge 
et  al.  [61],  Pita  and  Sundareasan  [62]  and  Hrenya  and  Sinclair  [63] 
incorporated  the  effects  of  gas  turbulence  into  their  models  by 
modifying  single-phase  turbulence  closures  to  account  for  the 
presence  of  particle  phase.  Samuelsberg  and  Hjertager  [64]  included 
the  gas  turbulence  by  using  LES  method  and  particle-particle 
interaction  using  the  kinetic  theory  approach.  Nieuwland  et  al. 
[65],  Balzer  et  al.  [66]  and  Neri  and  Gidaspow  [67]  conducted 
simulations  of  gas-solids  flow  in  the  circulating  fluidized  bed  using 
KTGF.  Armstrong  et  al.  [68]  applied  KTGF  to  study  the  wall-to-bed 
heat  transfer  in  a  bubbling  fluidized  bed.  Loha  et  al.  [69-71]  studied 
the  effect  of  gas-solid  drag  model,  the  solid  phase  wall  boundary 
condition  and  coefficient  of  restitution  on  the  hydrodynamic  beha¬ 
vior  of  a  bubbling  fluidized  bed  by  using  KTGF.  Papadikis  et  al. 
[72,73]  used  KTGF  to  study  the  effect  of  drag  model  and  particle  size 
in  fluidized  bed  system.  De  Wilde  and  Trujillo  [74]  applied  KTGF  to 
study  the  catalytic  cracking  of  gas  oil  in  a  rotating  fluidized  bed. 

In  KTGF,  the  kinetic  energy  associated  with  the  particle  velocity 
fluctuations  is  represented  by  a  pseudo-thermal  or  granular 
temperature  which  is  proportional  to  the  mean  square  of  the 
random  motion  of  particles.  A  separate  partial  differential  equa¬ 
tion  is  solved  for  the  granular  temperature  which  forms  the  basis 
for  the  turbulence  modeling  in  the  solid  phase.  Then,  the  solid 
phase  properties  like  solid  viscosity  and  pressure  are  described  as 
a  function  of  granular  temperature. 

The  granular  temperature  equation  is  represented  as 

=  (—PsI  +  Tg) 

:  V  v  s  +  V  ■  (k0s  V  0S) —  7 0s  “P  tpgs 

(47) 

where  (~psT +fg) :  V  v  s  is  the  generation  of  energy  by  solid  stress 
tensor.  V.(/cssV@s)  is  the  diffusion  of  energy  (kss  is  the  diffusion 
coefficient).  y0s  is  the  collisional  dissipation  of  energy  and  tpgs  is 
the  exchange  of  fluctuating  energy  between  phase. 

In  the  literature,  there  is  a  general  agreement  on  the  correlation 
of  solid  bulk  viscosity  and  solid  pressure  [56]  which  are  given 
below. 


2-(aspsSs )  4-  V  ■  (asps  V  S0S) 
dt 


where  p  is  the  coefficient  of  viscosity  which  depends  on  the 
thermodynamic  state  of  the  fluid  for  laminar  flow  and  have  an 
added  eddy  viscosity  for  turbulent  flow. 

Since  the  solid  phase  is  treated  as  continuous  fluid,  it  has 
similar  properties  to  a  fluid. 

Momentum  equation  for  the  solid  phase  is  written  as 

^(as/)s"v  S)  +  V  ■  (asps~Vs  V’s)=  —  asVp  —  Vps 4-  V  ■  7s4-asps~g  +Kgs(v  g  -  "v  s)+SgsV  s 

(45) 

where  ps  is  the  solid  pressure  and  ts  is  the  solid  phase  stress  tensor. 
The  stress  tensor  for  the  solid  phase  is  given  by 

ts  =  asps  rVvs  +  Vvs^  —  as  ^25  —  ^  (^  '  V  s)/  (46) 


—  3«sPsdsg0SS(l  +ess) 


1/2 


(48) 


Ps  —  asft0s  +  2/>s(l  +ess)a^goss0s  (49) 

where  ess  is  the  restitution  coefficient  and  goss  is  the  radial 
distribution  function. 

For  the  solids  shear  viscosity  different  descriptions  are  given  by 
different  authors  and  they  are  listed  in  Table  3.  For  example,  the 
equation  given  by  Gidaspow  [59]  does  not  account  for  the  inelastic 
nature  of  particles  in  the  kinetic  contribution  of  the  total  stress,  as 
Lun  et  al.  [56]  did.  The  solids  shear  viscosity  of  Syamlal  et  al.  [75] 
neglected  the  kinetic  or  streaming  contribution,  which  dominates  in 
dilute-phase  flow.  Flrenya  and  Sinclair  [63]  followed  Lun  et  al.  [56] 
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Table  3 

Solid  shear  viscosity. 

Equation 

Reference 

* [(£+*)  (i+(8/5y3r2)°,s°)] 

Lun  et  al.  [56] 

Bs  =  &pAg o(l  -eyf +^grr^[l  +§0  +e)(3e-  l)a,g0] 

Syamlal  et  al.  [75] 

,s  -  tfpMd  1  +  +*1 +e>*gb] 2 

Gidaspow  [59] 

_5v5 8  j  \f  1  1  ,  8a,  \  /'l+(8/5)l(37-2)a.g<A  .  768  2„  1 

Bs-  96  P^s[^  +  (Am,r/R),g,+  5  J{  2-„  j+ZSTOoJ 

Hrenya  and  Sinclair  [63] 

Table  4 

Gas-solid  drag  models. 


Equation 


cd  =  ^[1+°  15KRes)0687] 

CD=  (0.63  +  ^AL)  ,  Res  =  ^-J - fl 

\  s/R‘VA'r,  J  « 

vr,s  =  0.5^4-  0.06Res  +  \/(0.06Res)2  +  0.12(2B-A)+A2j , 

A  =  a414;  B  =  0.8a’ 28  for  af  <0.85;  B  =  a2  65  for  ag  >  0.85 
—  [w  +  ® V  1  ~  V  s|as“gr1'8 


Res  = 


dsag  Vs-  V  g \pg 


Kgs  —  [i^7+0-336j^|  v  s  -  v  g|as«g 


Res=- 


Kgs  =  1 50^|+1. 75^| | ,  ag<  0.8 
Rgs  =  I^D^I "^s  —  "^g  | «g  2  65 ,  «g  >  0.8 

[l  +  0.1 5(agRes)°  687] ,  Res  <  1 000 
0.44,  Res  >  1000 


Res  =- 


H 


Kgs  =  C(^+ 0.336) 

^ds|"Vs-"V  g\eK 

Res  = 


Kgs  =  1 50^+ 1. 75^| ts-5?g | ,  ag<  0.74 
%=iCD^|^s-^g|a> 

-°-5760+4<gg-o.r3i;+o.oo^  0.74  <  <%  <  0.82 


-0.0101 +- 


4(ag  —  0.7789r  +0.0040’ 
-31.8295+32.8295,  a,  >  0.97 


0.82  <  aE  <  0.97 


Reference 

Richardson  and  Zaki  [76] 
Wen  and  Yu  [77] 


Syamlal-O'Brine  [75] 


Gibilaro  et  al.  [82] 


Arastoopour  et  al.  [80] 


Gidaspow  [59] 


McKeen  and  Pugsley  [81] 


Yang  et  al.  [83] 


approach,  but  constrained  the  mean  free  path  of  the  particle  by  a 
dimension  characteristic  of  the  actual  physical  system. 

The  gas-solid  drag  model  plays  an  important  role  in  fluidized 
bed  modeling.  Different  gas-solid  drag  models  are  available  in  the 
literature  and  they  are  listed  in  Table  4.  An  earlier  drag  model  was 
proposed  by  Richardson  and  Zaki  [76].  Wen  and  Yu  [77]  proposed 
a  drag  model  by  extending  the  work  of  Richardson  and  Zaki  [76], 
The  model  proposed  by  Syamlal-O'Brien  [78]  which  was  based  on 
the  measurement  of  terminal  velocity  of  solid  particle  in  the 
fluidized  bed.  The  Gidaspow  mode  [59]  was  a  combination  of 
Wen  and  Yu  model  [77]  for  solid  volume  fraction  lower  than 
0.2  and  Ergun  model  [79]  for  solid  volume  fraction  equal  or  larger 


than  0.2.  Arastoopour  et  al.  [80]  proposed  a  drag  model  which 
gave  continuous  value  over  all  range  of  solid  volume  fraction. 
McKeen  and  Pugsley  [81]  proposed  a  drag  model  by  introducing  a 
constant  scale  factor  in  Gibilaro  et  al.  [82]  model  to  reduce  the 
drag  force  and  taking  into  account  the  inter  particle  cohesive  force 
on  particle  agglomeration.  To  investigate  the  dependence  of  drag 
coefficient  on  structure  parameters,  Yang  et  al.  [83]  adopted  a 
structure-dependent  drag  model  based  on  the  energy  minimiza¬ 
tion  multi-scale  (EMMS)  approach. 

The  Euler-Euler  two-fluid  model  in  combination  with  KTGF  is 
widely  applied  for  investigating  the  hydrodynamics  of  fluidized 
bed  system.  But,  it  is  difficult  to  couple  the  complicated  chemical 
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reactions  with  the  fluidized  bed  hydrodynamics  due  to  the  complex 
mechanism  of  heat  transfer  and  chemical  reaction  need  to  be 
modeled  in  fluidized  bed  gasification  process.  To  model  the  chemical 
reactions,  more  careful  consideration  has  to  be  given  on  the  solution 
of  large  numbers  of  energy  and  species  transport  equations  and 
nonlinear  source  terms  of  complicated  chemical  reactions.  For 
modeling  the  fluidized  bed  system  with  chemical  reaction  by  using 
Euler-Euler  model,  in  addition  to  the  continuity  and  momentum 
equation,  species  conservation  equation  for  each  gas  species  and 
energy  equation  for  gas  and  solid  phase  is  required  to  solve. 

The  species  transport  equation  is  given  by 

—(PgdgYg  j)  +  V(pgCig  v  gYgj)  =  —  V  ■  (otg]g  i)  +  agRg,i  +  Rs,i  (50) 

where  Jgi,  Rgi  and  Rsi  are  the  diffusion  flux  of  species  i  in  gas 
phase,  the  net  rate  of  production  of  homogeneous  species  i  and  the 
homogeneous  reaction  rate,  respectively. 

In  species  transport  equations  of  gas  phase,  mass  diffusion 
coefficient  are  used  to  calculate  the  diffusion  flux  of  chemical 
species  in  turbulent  flow  using  modified  Fields  law 


(51) 


where  aY  is  the  Schmidt  number,  which  is  set  as  0.7.  The  diffusion 
coefficient  of  the  mixture,  Dm  i  is  calculated  from  the  binary  mass 
diffusion  coefficient  D, ,j  as  follows: 


Dm,i  = 


1-X, 
I,  •-  iXj 


(52) 


The  energy  transport  equations  are  solved  for  the  specific  enthalpy 
of  gas  phase  and  solid  phase,  which  take  the  form: 

i)t-(agPgHgl  +  V  ■  (agpg  v  gHg)  =  V(kgTg)+QgS  +  SgsHs  (53) 

and 

— (^sPsFfs)  + V  ■  (atsPs  v  SHS)  =  V(/csTs)  +  QSg  +  SSgFfs  (54) 


where  H,  k,  Q  are  the  specific  enthalpy,  the  gas  mixture  thermal 
conductivity  and  the  intensity  of  heat  exchange  between  the  gas  and 
solid  phases,  respectively.  The  third  term  on  the  right  hand  side  is  the 
heat  transfer  in  that  the  solid  phase  changed  into  gas  phase. 

The  specific  enthalpy  is  defined  by 


H=  £  y,H, 
!'=  1 


(55) 


where  Hi  is  the  enthalpy  for  each  chemical  species  in  the  mixture 
and  considers  both,  thermal  and  chemical  enthalpy 

H(=  fT  CPidT+AHfi  (56) 

Jt0 

where  T0,  CPi  and  A Hfi  are  the  reference  temperature,  the  heat 
capacity  at  constant  pressure  for  the  ith  species,  and  the  enthalpy 
of  formation  for  the  ith  species  in  the  standard  state. 

The  above  conservation  equations  are  solved  by  using  CFD 
technique  to  get  the  gas  and  solid  velocity  distribution,  volume 
fraction  distribution,  gas  composition  distribution  and  tempera¬ 
ture  distribution  inside  the  gasifier  and  their  variation  with  time. 
A  typical  solution  procedure  for  Euler-Euler  two-fluid  modeling  of 
fluidized  bed  gasification  process  using  Phase  Couple  SIMPLE 
algorithm  is  shown  in  Fig.  8. 

Very  few  studies  on  fluidized  bed  gasification  process  using 
Euler-Euler  model  are  available  in  the  literature.  Dou  et  al.  [84] 
and  Dou  and  Song  [85]  studied  the  hydrogen  production  from 
steam  reforming  of  glycerol  in  a  lab-scale  fluidized  bed  reactor  by 
using  the  Euler-Euler  model.  The  simulation  was  performed  for  a 
two-dimensional  computational  domain.  The  simulated  flow 


Fig.  8.  Flowchart  for  solution  procedure  of  Euler-Euler  CFD  model. 


Table  5 

Comparison  between  simulated  and  experimental  data  for  steam  reforming  of 
glycerol  in  fluidized  bed  reactor  [87], 


Experimental  data  Simulation 

Inlet  gas  velocity  (m/s) 

0.5 

0.5 

Steam  to  carbon  molar  ratio  (S/C) 

2:1 

2:1 

Catalyst 

NiO/Al203 

NiO/Al203 

Temperature  (  C) 

600 

600 

H2  (vol%) 

59.1 

52.6 

C02  (vol%) 

27.2 

33.5 

CH4  (vol%) 

5.9 

13.9 

CO  (vol%) 

7.8 

0.0  (3.5  x  10~15) 

distribution  exhibited  a  more  heterogeneous  core-annuls  structure 
of  gas-solid  flow  which  led  to  back  mixing  and  internal  circula¬ 
tion.  This  flow  structure  suggested  that  the  bed  should  be  agitated 
to  maintain  satisfactory  fluidizing  conditions.  The  glycerol  conver¬ 
sion  and  H2  production  variation  with  inlet  gas  velocity  was 
studied.  The  simulated  results  were  also  validated  with  experi¬ 
mental  results  as  shown  in  Table  5.  The  difference  between 
simulated  and  experimental  gas  composition  was  less  than  8%. 
But,  the  model  predicted  very  less  (3.5  x  10~15%)  amount  of  CO. 
Though,  it  was  mentioned  that  the  error  incurred  in  calculation 
derived  from  several  sources  such  as  flow  and  reaction  rate  model, 
numerical  treatment  of  initial  and  boundary  conditions  and  other, 
but,  the  specific  reason  for  predicting  very  less  amount  of  CO 
composition  was  not  described.  Papadikis  et  al.  [86,87]  studied  the 
fast  pyrolysis  of  biomass  in  a  lab-scale  bubbling  fluidized  bed 
reactor  by  using  Euler-Euler  model.  The  Eulerian  approach  was 
used  to  model  the  bubbling  behavior  of  sand.  The  discrete 
description  was  used  for  modeling  the  biomass  particle  injected 
inside  the  reactor  and  the  motion  of  the  particle  was  computed 
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using  local  volume  fraction  dependent  drag  laws.  2-D  and  3-D 
simulations  were  run.  It  was  showed  that  in  the  3-D  simulation  a 
bubble  was  formed  close  to  the  feeding  point  of  the  reactor  and 
the  drag  force  was  quite  reduced  at  the  injection  time  compared 
with  the  2-D  case.  The  simulated  hydrodynamic  behavior  of  the 
bed  at  different  time  where  biomass  particle  position  was  indi¬ 
cated  by  black  spot  is  shown  in  Fig.  9.  The  model  gave  important 
information  like  hydrodynamic  of  bed  with  biomass  particle 
position,  degradation  and  evolution  of  tar  with  time  but  no 
experimental  validation  was  provided.  Xue  et  al.  [88]  also  simu¬ 
lated  the  fast  pyrolysis  of  biomass  in  a  lab-scale  fluidized  bed 
reactor  by  using  Euler-Euler  model.  Here,  gas-particle  flow  was 
modeled  with  a  multi-fluid  description  (gas,  sand  and  biomass). 
Results  indicated  that  biomass  particle  size  and  superficial  gas 
velocity  influenced  the  tar  yield  and  residence  time  considerably 
with  a  fixed  bed  height.  The  simulated  yields  of  bio-oil,  char  and 
non-condensable  gas  were  compared  with  the  experimental  data 
for  pure  cellulose  and  red  oak  as  shown  in  Table  6.  The  predicted 
results  showed  good  quantitative  agreement  with  the  experimen¬ 
tal  data  for  the  yield  of  each  product.  But,  the  model  over¬ 
predicted  the  unreacted  biomass  (residue)  particularly  for  smaller 
particles.  Therefore,  it  was  suggests  that  a  particle  size  distribution 
may  need  to  be  considered  to  improve  the  simulation.  Yu  et  al.  [89] 
simulated  the  lab-scale  bubbling  fluidized  bed  coal  gasifier  by 
using  Euler-Euler  model.  Their  model  was  two  dimensional  and 
the  gasification  kinetic  was  described  by  15  species  and  11 
chemical  reactions  which  include  pyrolysis,  homogeneous  reac¬ 
tions  and  heterogeneous  reactions.  The  devolatalization  and  dry¬ 
ing  was  assumed  as  instantaneous  process  in  the  feed  zone. 
Pressure,  temperature,  velocity,  volume  fraction  of  gas  composi¬ 
tion  were  predicted.  The  simulated  distribution  of  gas  composi¬ 
tions  within  the  gasifier  was  plotted  for  six  different  operating 
conditions  and  two  of  them  are  shown  here  in  Fig.  10.  It  was 
observed  that  overall  trend  of  each  composition  was  in  consistence 
with  the  operating  conditions.  For  example,  the  concentrations  of 
C02  and  CO  increased  along  the  height  of  the  gasifier,  while  H2  and 
CH4  go  up  at  first  and  then  drop  down  to  the  top  of  the  reactor  due 
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Fig.  9.  Fluidized  bed  hydrodynamics  with  biomass  particle  position  [89], 


to  devolatilization  and  equilibrium  of  water-gas-shift  reactions. 
This  showed  that  the  Euler-Euler  model  gave  more  exact  predic¬ 
tion  which  could  not  be  described  by  two-phase  flow  model.  For 
validation,  the  predicted  outlet  gas  compositions  for  six  different 
operating  conditions  were  compared  with  the  experimental  result 
from  the  literature  as  shown  in  Fig.  11.  It  was  observed  that  the 
predicted  results  agreed  well  to  the  experiment  with  5%  error  in 
C02  and  N2  molar  fraction  and  within  20  %  error  for  other  results. 
Wang  et  al.  [90]  extended  the  above  simulation  in  three  dimen¬ 
sional  geometry  to  take  into  account  the  influences  of  bed  walls 
and  geometry  on  the  flow  pattern  and  other  aspects  in  3D. 
Armstrong  et  al.  [91]  also  simulated  the  same  lab-scale  bubbling 
fluidized  bed  coal  gasifier  using  Euler-Euler  model.  They  have 
introduced  the  limestone  calcinations  to  their  model.  Their  model 
considered  two  solid  phases  one  for  limestone  and  other  for  char 
which  enabled  to  simulate  the  segregation  of  lower  density 
particles  (char)  to  the  top  of  the  bed.  The  hydrodynamic  investiga¬ 
tion  revealed  the  bubble  formation  exogenously  and  endogenously 
as  a  result  of  the  reaction  kinetics  as  shown  in  Fig.  12.  A  detailed 
analysis  of  gasification  process  was  presented  from  the  distribu¬ 
tion  of  gas  composition,  temperature  and  reaction  rate  inside  the 
gasifier.  The  model  results  considering  inert  limestone  and  lime¬ 
stone  calcinations  were  compared  to  the  experimental  data  from 
the  literature  as  shown  in  Table  7.  The  model  gave  reasonably  good 
representation  of  experimental  compositions  as  a  whole.  However, 
unlike  experiment,  model  1  predicted  higher  CO  than  C02.  They 
expected  that  this  was  due  to  fact  that  experimental  bed  consisted 
of  limestone  whereas  the  model  was  formed  of  both  char  and 
limestone.  However,  model  2  gave  higher  C02  compared  to  model 
1  due  to  the  lower  initial  temperature  of  the  bed.  The  gas 
composition  did  not  change  significantly  by  considering  limestone 
calcinations  due  to  slow  conversion  rate  of  limestone  decomposi¬ 
tion.  The  above  simulation  was  also  extended  to  determine  the 
effects  of  different  bed  ratios  and  over  a  longer  period  of  time  [92], 
The  above  literature  showed  that,  the  Euler-Euler  model  is  an 
efficient  way  to  study  the  fluidized  bed  gasification  process  where 
the  gas  and  the  solid  phases  are  treated  as  the  interpenetrating 
continuum  and  detail  momentum  equations  are  solved  for  both 
the  phases  in  combination  with  reaction  kinetics  to  get  a  detailed 
insight  into  the  gasifier.  In  contrast  to  the  two-phase  flow  model, 
the  Euler-Euler  model  can  simulate  the  dense  bottom  bed  and 
freeboard  region  both  simultaneously  and  the  three  dimensional 
flows  can  be  solved.  It  can  provide  much  more  detailed  informa¬ 
tion  compared  to  the  two-phase  flow  model  which  are  very  much 
useful  for  design  and  optimization.  The  Euler-Euler  model  can 
predict  the  bubble  formation,  bubble  rice  through  the  bed,  inter¬ 
action  with  other  bubbles  and  growth  in  size  and  eruption  at  the 
upper  surface  of  the  bed  during  fluidized  bed  gasification  process. 
It  can  provide  the  axial  and  radial  variation  of  gas  and  solid  volume 
fractions,  gas  composition,  temperature  and  reaction  rates  inside 
the  gasifier.  In  addition  to  the  effect  of  operating  parameters,  other 
effects  like  position  of  fuel  inlet,  gas  inlet,  outlet,  geometry  of  the 


Table  6 

Comparison  between  simulated  and  experimental  product  yield  for  fluidized  bed  fast  pyrolysis  of  biomass  [88], 


Method 

Bio-oil 

Char 

Non-considerable  gas 

Residual 

Temperature  (  C) 

Product  yields  (wt%)  of  pure  cellulose  pyrolysis  from  experiment  and  simulation 

Experiment 

82.1 

2.2 

12.4 

- 

500 

Simulation 

82.2 

3.3 

13.9 

0.5 

493 

Product  yields  (wt%)  of  a  type  of  red  oak  pyrolysis  from  experiment  and  simulation 

Experiment 

71.7  ±  1.4 

13.0  +  1.5 

20.5  + 1.3 

- 

500 

Simulation  (dp= 250  pm,  2D) 

60.5 

12.3 

16.2 

11.0 

497 

Simulation  (dp= 325  pm,  2D) 

62.4 

14.1 

17.3 

6.2 

498 

Simulation  (dp= 325  pm,  2D) 

61.5 

12.9 

16.5 

8.7 

499 

Simulation  (dp=400  pm,  2D) 

63.4 

15.1 

18.1 

3.4 

499 
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Case  I 


Fig.  10.  Concentration  distribution  of  H2,  C02,  N2,  CH4  and  CO  for  different  operating  conditions:  Case  1 :  8.0  kg/h  coal  feed,  21.9  kg/h  air  supply,  4.6  kg/h  steam  supply,  420  C 
air  and  steam  temperature,  855  “C  reactor  temperature.  Case  2:  8.0  kg/h  coal  feed,  17.0  kg/h  air  supply,  4.6  kg/h  steam  supply,  413  "C  air  and  steam  temperature,  and  812  C 
reactor  temperature  [84], 


gasifier,  initial  bed  height  etc.  can  also  be  studied.  In  Euler-Euler 
model,  the  mixing  and  segregation  between  inert  bed  material 
and  solid  fuel  particle  inside  the  gasifier  can  also  be  studied  by 
taking  two  different  solid  phases  with  different  material  proper¬ 
ties.  The  Euler-Euler  model  is  computationally  less  expensive 
compare  to  the  Euler-Lagrange  model  because  it  assumes  particles 
as  continuum  instead  of  tracking  each  and  individual  particle  and 
hence  can  be  used  for  modeling  large  scale  fluidized  bed  gasifier. 
But,  due  to  the  approximation  of  particle  flows  as  fluid,  success  of 
the  Euler-Euler  model  depends  on  the  proper  description  of 
particle  viscosity  and  pressure.  Besides,  to  simulate  particle  size 
distribution,  more  than  one  pseudo-fluid  has  to  be  included  and 
all  conservation  equations  are  needed  to  be  solved  for  each 
pseudofluid  which  represents  one  particle  size.  Therefore,  gener¬ 
ally  Euler-Euler  modeling  is  performed  for  one  or  two  particle  size 
not  for  a  wide  range  of  particle  size  distribution.  The  Euler-Euler 
model  cannot  provide  information  about  the  residence  time  of 
individual  reacting  particles.  Another  drawback  of  Euler-Euler 
model  is  that,  the  particle  degradation  and  shrinkage  during 


gasification  process  cannot  be  modeled  here.  Therefore,  to  get 
more  particle  level  information  of  fluidized  bed  gasification 
process,  Euler-Lagrange  modeling  is  required. 

4.3.  Euler-Lagrange  model 

In  Euler-Lagrange  model,  the  gas  phase  is  treated  as  a  con¬ 
tinuum  and  time-averaged  Navier-Stoke  equations  are  solved  for 
the  gas  phase,  while  the  solid  phase  is  solved  by  tracking  each  and 
individual  particle  in  Lagrangian  frame  of  reference.  The  solid 
phase  exchange  mass,  momentum  and  energy  with  the  gas  phase. 
In  Euler-Lagrange  model,  for  particle  description,  discrete  element 
model  (DEM)/  discrete  particle  model  (DPM)  is  commonly  applied. 
In  DEM/DPM  model,  the  collision  between  particles  may  be  either 
based  on  soft-sphere  approach  or  hard-sphere  approach. 

In  a  hard-sphere  approach,  the  trajectories  of  the  particles  are 
determined  by  momentum-conserving  binary  collisions.  The  interac¬ 
tions  between  particles  are  assumed  to  be  pair-wise  additive  and 
instantaneous.  In  the  simulation,  the  collisions  are  processed  one  by 


Molar  traction  (%)  Molar  fraction  (%)  Molar  fraction  (%) 
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Fig.  11.  Comparison  between  predicted  and  experimental  outlet  gas  composition  for  different  cases  [84], 
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Fig.  12.  Volume  fraction  distribution  within  the  bed  for  (A)  gases,  (B)  limestone  and  (C)  char  [90]. 
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Table  7 

Comparison  of  model  results  with  inert  limestone  and  calcinating  limestone  and  experimental  result  [91], 


CO 

co2 

H2 

ch4 

n2 

h2o 

Tar 

Model  1 

Inert  limestone 

0.13464 

0.09150 

0.06226 

0.00020 

0.54675 

0.16428 

0.00037 

Limestone  calcinations 

0.13249 

0.09279 

0.06094 

0.00020 

0.54671 

0.16649 

0.00037 

Difference 

-0.00214 

0.00130 

-0.00132 

0.00000 

-0.00003 

0.00220 

0.00000 

Model  2 

Inert  llimestone 

0.10135 

0.09742 

0.05230 

0.00014 

0.52268 

0.22584 

0.00026 

Limestone  calcinations 

0.10268 

0.09818 

0.05422 

0.00014 

0.52379 

0.22073 

0.00026 

Difference 

0.00132 

0.00076 

0.00192 

0.00000 

0.00111 

-0.00511 

0.00000 

Experiment 

Exp.  1 

0.1094 

0.1931 

0.0853 

0.0084 

0.6037 

N/A 

N/A 

Exp.  2 

0.1059 

0.1838 

0.0884 

0.0107 

0.6110 

N/A 

N/A 

one  according  to  the  order  in  which  the  events  occur.  For  not  too 
dense  systems,  the  hard-sphere  models  are  considerably  faster  than 
the  soft-sphere  models.  Occurrence  of  multiple  collisions  at  the  same 
instant  cannot  be  taken  into  account  in  hard  sphere  model  [93], 
Campbell  and  Brennen  [94]  reported  the  first  hard-sphere  simulation 
used  to  study  granular  systems.  Since  then,  the  hard-sphere  models 
are  applied  to  study  a  wide  range  of  complex  granular  systems. 
Hoomans  et  al.  [95]  used  the  hard-sphere  model  in  combination  with 
a  CFD  approach  for  the  gas-phase  conservation  equations  to  study 
gas-solid  two  fluidized  beds.  By  using  this  model,  they  studied  the 
effect  of  particle-particle  interaction  on  bubble  formation  and  the 
particle  segregation  induced  by  particle  size  differences  and  density 
differences  [96],  This  model  has  been  further  used  to  study  high- 
pressure  fluidization  [97,98],  circulating  fluidized  [99],  spout-fluid 
beds  [100,101]  and  particle  flows  through  contractions  [102].  Similar 
simulations  are  also  carried  out  by  many  other  research  groups.  For 
example,  Dahl  et  al.  [103]  and  Dahl  and  Hrenya  [104,105]  applied  a 
hard-sphere  model  to  study  segregation  in  continuous  size  distribu¬ 
tions.  Ouyang  and  Li  [106,107]  developed  a  slightly  different  version  of 
this  model.  Helland  et  al.  [108]  developed  a  DPM  in  which  hard-sphere 
collisions  are  assumed,  but  a  time-driven  scheme  is  used  to  locate  the 
collisional  particle  pair.  Effect  of  the  gas  turbulence  has  also  been  taken 
into  account  in  some  hard-sphere  models  [109-114].  At  high  particle 
number  densities  or  low  coefficient  of  normal  restitution,  the  collision 
will  lead  to  a  dramatical  decrease  in  kinetic  energy.  This  is  the  so  called 
inelastic  collapse,  in  which  regime  the  collision  frequencies  diverge  as 
relative  velocities  vanish.  In  that  case,  the  hard  sphere  method 
becomes  useless.  In  more  complex  situations,  the  particles  may 
interact  via  short  or  long-range  forces  and  the  trajectories  are 
determined  by  integrating  Newtonian  equation  motion.  The  soft- 
sphere  method  originally  developed  by  Cundall  and  Strack  [115]  was 
the  first  granular  dynamics  simulation  technique  published  in  the 
open  literature.  Soft-sphere  models  use  a  fixed  time  step  and  conse¬ 
quently  the  particles  are  allowed  to  overlap  slightly.  The  contact  forces 
are  subsequently  calculated  from  the  deformation  history  of  the 
contact  using  a  contact  force  scheme.  The  soft-sphere  models  allow 
for  multiple  particle  overlap  although  the  net  contact  force  is  obtained 
from  the  addition  of  all  pair-wise  interactions.  The  soft-sphere  models 
are  essentially  time  driven,  where  the  time  step  should  be  carefully 
chosen  in  the  calculation  of  the  contact  forces.  The  soft-sphere  models 
that  can  be  found  in  literature  mainly  differ  from  each  other  with 
respect  to  the  contact  force  scheme  that  is  used.  A  review  of  various 
popular  schemes  for  repulsive  inter-particle  forces  was  presented  by 
Schafer  et  al.  [116].  Walton  and  Braun  [117]  developed  a  force  model 
which  uses  two  different  spring  constants  to  model  the  energy 
dissipation  in  the  normal  and  tangential  directions  respectively.  In 
the  force  scheme  proposed  by  Langston  et  al.  [118],  a  continuous 
potential  of  an  exponential  form  was  used,  which  contains  two 
unknown  parameters:  the  stiffness  of  the  interaction  and  an  interac¬ 
tion  constant.  A  two-dimensional  soft-sphere  approach  was  first 


applied  to  gas-fluidized  beds  by  Tsuji  et  al.  [119],  where  the  linear 
spring/dashpot  model  similar  to  the  one  presented  by  Cundall  and 
Strack  [115]  was  employed.  Kawaguchi  et  al.  [120]  extended  this 
model  to  three  dimensions  as  far  as  the  motion  of  the  particles  is 
concerned.  Yu  and  co-workers  [121-123]  independently  developed  a 
two-dimensional  model  of  a  gas-fluidized  bed.  Flowever,  in  their 
simulations  a  collision  detection  algorithm  that  is  normally  found  in 
hard-sphere  simulations  was  used  to  determine  the  first  instant  of 
contact  precisely.  Feng  and  Yu  [124]  and  Feng  et  al.  [125]  applied  this 
model  to  study  segregation  processes  of  a  binary  mixture.  Iwadate  and 
Horio  [126]  and  Mikami  et  al.  [127]  developed  a  model  based  on  Tsuji 
et  al.  [119]  where  they  incorporated  Van  der  Waals  forces  to  simulate 
fluidization  of  cohesive  particles.  Further  studies  of  the  influence  of  gas 
and  particle  properties  for  Geldart  A  particles  were  performed  by  Ye 
et  al.  [128,129]  and  Pandit  et  al.  [130].  Kafui  et  al.  [131]  developed  a 
DPM  based  on  the  theory  of  contact  mechanics,  thereby  enabling  the 
collision  of  the  particles  to  be  directly  specified  in  terms  of  material 
properties  such  as  friction,  elasticity,  elastoplasticity  and  auto¬ 
adhesion.  Limtrakul  et  al.  [132]  used  the  soft-sphere  model  with  mass 
transfer  and  chemical  reactions  to  study  the  decomposition  of  ozone 
on  catalyst  coated  particles  in  a  two-dimensional  fluidized  bed. 
Kuwagi  et  al.  [133]  coupled  the  soft-sphere  model  with  a  model  for 
the  description  of  metallic  solid  bridging  by  surface  diffusion  mechan¬ 
isms  including  the  effect  of  surface  roughness.  Oevermann  et  al.  [134] 
studied  Euler-Lagrange/DEM  simulation  using  soft-sphere  method  to 
study  wood  gasification  in  bubbling  fluidized  bed.  The  soft-sphere 
DEM  model  combining  with  CFD  was  also  used  in  the  literature 
[135,136]  to  study  the  heat  transfer  in  bubbling  fluidized  beds. 
Recently,  In  the  following  section,  only  soft-sphere  DEM/DPM  is 
discussed  because  it  can  in  principle  handle  any  situation  (dense 
regimes,  multiple  contacts)  and  also  additional  interaction  force. 
Mathematical  equations  used  for  describing  the  soft-sphere  DEM/ 
DPM  are  given  below. 

For  each  particle,  the  linear  momentum  equations  is  given  by 


dVj 

'dt 


d  n 
1 dt 2 


=  Fr, 


i  +  Fe] 


(57) 


where  the  RF1S  is  the  total  force  on  particle,  which  has  three  basic 
contributions: 


(i)  The  total  contact  force  Fcontacti,  which  is  the  sum  of  individual 
contact  forces  exerted  by  all  other  particles  being  in  contact 
with  the  particle  i,  which  are  divided  into  a  normal  and  a 
tangential  component: 

F contact,!  =  Fcontact,N,i  Fcon[aC[  y  i  (58) 


(ii)  The  total  external  force  Fext 
Fext.i  =  Fgi  4-  Fdj  +  Fpi 


(59) 
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where  Fgi  is  the  gravitational  force,  Fdi  is  the  drag  force 
exerted  by  the  surrounding  gas  phase  and  Fpi  is  the  pressure 
force. 

(iii)  The  sum  of  all  other  particle-particle  forces  Fpp,a,  which  can 
include  short  range  cohesive  force  form  the  van  der  Walls 
interaction  between  the  molecules,  as  well  as  long  range 
electrostatic  forces. 

The  angular  momentum  equation  of  particle  is  given  by 


where  /,•  is  the  moment  of  inertia  of  particle  i  and  7j-  is  the 
torque  which  depends  only  on  the  tangential  component  of 
the  individual  contact  forces. 

The  heat  balance  for  individual  particle  is  given  by 

ITliCpj-^^^  Qgp  Qpp  ~F  Qrad  ~F  Qr  (61) 

where  Qgp  is  the  gas-solid  convection  heat  transfer,  Qpp  is  the 
particle-particle  conduction  heat  transfer,  Qrad  is  the  radiation 
heat  transfer  and  QR  is  the  heat  transfer  due  to  heterogeneous 
chemical  reactions. 

The  gas-phase  mass  and  momentum  equations  are  as  follows: 
^(agpg)  +  V  ■(agpg~vg)  =  Sg  (62) 

—(agPg  V  g)  +  V  ■  (agpg  V  g  V  g)  =  — ctgVPg  +  V  ■  (agTg)+agpg  g  +  Sp 

(63) 

where  a  is  the  volume  fraction,  p  is  the  density,  "v"  is  the  velocity 
vector  of  gas  phase.  Sg  is  the  mass  source  term  due  to  chemical 
reactions.  zg  is  the  viscous  stress  tensor  and  Sp  is  the  momentum 
source  due  to  inter-phase  interaction.  For  gas-solid  phase  in 
fluidized  bed,  two-way  coupling  is  required.  The  Sp  is  computed 
by  adding  up  the  drag  force  of  the  particles  located  in  fluid  cell 
which  is  given  by 

1  np 

Sp  =  ts —  2  Fd  (64) 

v  cell  i  =  1 

where  V ced  is  the  volume  of  the  fluid  cell. 

The  gas-phase  energy  equation  is  given  by 

— (agPgCpgTg)-l- V  ■  (agpg  v  gCpgT)=  V  ■  (agkgVTg)-t-SQcv  +  SQR  +  Sf, 

(65) 

where  Tg  is  the  temperature,  6pg  is  the  heat  capacity,  kg  is  the 
thermal  conductivity  of  the  gas  phase.  Sq  Cv  is  the  heat  source  due 
to  gas-solid  convective  heat  transfer,  Sq  R  is  the  heat  source  due  to 
chemical  reactions  and  Sh  is  the  heat  transported  by  the  mass 
source  Sg. 

Euler-Lagrange  DEM/DPM  allows  solution  for  flows  with  a 
wide  range  of  particle  types,  sizes,  shapes  and  velocities.  It  is  able 
to  look  at  individual  discrete  particles  and  their  interaction  with  its 
local  environment.  But,  due  to  high  collision  frequency  for  volume 
fractions  above  5%  and  the  computational  complexity  of  calculat¬ 
ing  dense  particle-particle  interactions,  DPM  calculations  have 
been  limited  to  on  the  system  having  comparatively  less  numbers 
of  particles  [137],  Generally,  Euler-Lagrange  DEM/DPM  is  used  to 
study  the  combustion  and  gasification  of  fuel  in  freeboard  area 
where  the  particle  volume  fraction  is  less.  There  has  been  limited 
study  on  the  simulation  of  dense  gas-solid  flow  coupling  with 
chemical  reactions  by  using  Euler-Lagrange  DEM/DPM.  Silaen  and 
Wang  [138]  and  Watanabe  and  Otaka  [139]  applied  Euler- 
Lagrange  DEM/DPM  in  entrained  flow  coal  gasifier.  Grabner  et  al. 
[140]  used  Euler-Lagrange  DEM/DPM  at  circulating  fluidized  bed 


conditions  but  did  not  take  particle  collisions  into  account. 
Recently,  Bruchmuller  et  al.  [141]  studied  the  thermochemical 
degradation  of  biomass  inside  a  lab-scale  (38.1  mm  diameter  and 
460  mm  height)  fast  pyrolysis  fluidized  bed  reactor  by  using 
Euler-Lagrange  DEM/DPM.  They  have  tacked  0.8  million  individual 
discrete  sand  and  biomass  particles.  This  made  it  possible  to  look 
at  accurate  and  detailed  multi-scale  information  (i.e.,  any  desired 
particle  property,  trajectory,  and  particle  interaction)  over  the 
entire  particle  life  time  which  revealed  that  the  overall  thermo¬ 
chemical  degradation  process  of  biomass  was  influenced  by  local 
flow  and  particle  properties.  Multiprocessor  simulation  with  high- 
end  computer  (2.3  GHz,  192  cores  in  total  distributed  on  48  quad 
cores)  was  carried  out  which  took  320,000  CPU  hours  to  run  5  s 
real  time. 

The  Euler-Lagrange  DEM/DPM  is  the  most  comprehensive 
model  developed  till  date  to  study  the  fluidized  bed  system. 
Because,  it  has  the  highest  potential  to  realistically  reveal  thermo¬ 
chemical  processes  in  granular  multiphase  flow  applications. 
DEM/DPM  models  are  able  to  accurately  predict  interparticle  and 
interphase  exchange  of  mass,  momentum  and  energy  at  the 
individual  particle  level  without  any  approximation.  Therefore, 
trajectories,  temperature,  composition  and  any  other  additional 
particle  information  are  more  reliably  and  more  naturally  included 
which  are  essential  to  study  effects  like  thermochemical  degrada¬ 
tion,  shrinkage,  breakage,  segregation,  mixing  and  entrainment. 
But,  the  computational  time  and  cost  involved  in  Euler-Lagrange 
DEM/DPM  simulation,  restricting  the  applicability  of  the  model 
particularly  for  simulating  large-scale  fluidized  bed  reactor  which 
involves  more  particles  and  reaction  kinetic.  Besides,  due  to  its 
complexity,  Euler-Lagrange  DEM/DPM  simulations  are  mostly 
performed  in  2-D  or  quasi-3-D. 

Very  recently,  Euler-Lagrange  Computational  Particle  Fluid 
Dynamics  (CPFD)  model  is  employed  for  calculating  dense  particle 
flow.  The  CPFD  numerical  methodology  incorporates  the  multi- 
phase-particle-in-cell  (MP-PIC)  method  [142,143],  where  the  fluid 
phase  is  solved  by  using  an  Eulerian  computational  grid  and  the 
solids  are  modeled  using  Lagrangian  computational  particles.  Here 
a  computational  or  numerical  particle  is  accounted  for  via  an 
ensemble  of  particles  displaying  the  same  properties  such  as 
chemical  composition,  size  and  density.  Hence,  the  number  of 
particles  needed  to  solve  can  be  reduced  from  billions  to  millions 
and  this  model  can  be  applied  to  simulate  lean  as  well  as  dense 
particle  phase  flow.  Therefore,  this  model  is  in  between  Euler- 
Euler  two-fluid  model  and  Euler-Lagrange  DEM/DPM.  Andrews 
and  O'Rourke  [142]  first  proposed  the  MP-PIC  method  and 
demonstrated  the  method  with  one-dimensional  simulations 
compared  to  analytical  solutions  and  experimental  data.  In  CPFD 
three-dimensional  forces  on  each  particle  is  considered  which 
include  fluid  drag,  gravity,  static-dynamic  friction,  particle  colli¬ 
sion  and  possibly  other  forces.  Using  this  method,  the  particle 
stress  gradient,  which  is  difficult  to  calculate  for  each  particle  in 
dense  flow,  is  calculated  as  a  gradient  to  the  grid,  fully  coupled  to 
the  other  particle  and  gas  acceleration  terms  and  is  then  inter¬ 
polated  to  the  discrete  particles  [144].  In  CPFD,  cell  averaged 
chemistry  is  used  to  calculate  the  reaction  kinetics.  Average 
properties  for  the  particle  phase  in  the  chemical  rate  equations 
are  calculated  by  interpolating  discrete  computational  particle 
properties  to  the  grid.  The  reaction  rates  are  calculated  in  each 
grid  cell  by  solving  a  set  of  ordinary  differential  equation.  The 
mathematical  equations  used  for  describing  the  fluid  and  solid 
phase  in  CPFD  method  are  given  below. 

The  volume  averaged  fluid  phase  mass  and  momentum  equa¬ 
tions  are  given  by 

d  _ >.  r 

-^afPf)  +  v  '  (afPf  v  f)=Pfs  (66) 
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dt 


(afPf  V  f)  +  W  ■  (ccfPf  VfVf) 


-  Vp-F+afPfg  +  V  ■  afrf 


(67) 


where  Vf  is  the  fluid  velocity  and  aj  is  the  fluid  volume  fraction, 
Pf  is  the  fluid  density,  p  is  the  fluid  pressure,  tf  is  the  fluid  stress 
tensor  and  g  is  the  gravitational  acceleration.  F  is  the  rate  of 
momentum  exchange  per  unit  volume  between  the  fluid  and  the 
solid  phase.  The  fluid  mass  production  rate  per  unit  volume  from 
solid-fluid  chemistry  is  pjs. 

The  constitutive  equation  for  the  fluid  stress  r/  is  given  by 


rf,ij  =  F 


/dVi 

\dXj  dXi) 


2  dV[ 
"3 


(68) 


where  p  is  the  coefficient  of  viscosity  which  depends  on  the 
thermodynamic  state  of  the  fluid  for  laminar  flow  and  have  an 
added  eddy  viscosity  for  turbulent  flow  based  on  Smagorinsky 
turbulence  model  [145], 

The  turbulence  viscosity  is  given  by 


Ft  —  C2PfA2 


(69) 


where  C  is  the  subgrid  scale  eddy  coefficient  and  A  is  the 
length  scale. 

The  fluid  phase  species  equation  is  given  by 
tjt-(afPfYf,i)  +  V  ■  ( ajpfYfjV’f )  =  V  ■  ( PfafDS/Yfj)  +  Smichem  (70) 


where  Vy  ,■  is  the  mass  fraction  of  each  fluid  species  and  Smi  chem  is 
the  net  production  rate  of  species  i  due  to  fluid  phase  chemical 
reactions.  The  coefficient  D  represents  the  turbulent  mass  diffu- 
sivity  which  is  related  to  the  viscosity  by  the  Schmidt  number 
[p/PfD  =  Sc). 

The  fluid  phase  energy  equation  in  terms  of  the  enthalpy  is 
given  as 


where,  R  is  the  universal  gas  constant,  Tf  is  gas  mixture  tempera¬ 
ture  and  Mi  is  the  molecular  weight  of  gas  species  i. 

The  dynamics  of  particle  phase  is  calculated  by  solving  a 
transport  equation  for  the  particle  distribution  function  (PDF)  /. 
It  is  assumed  that  /  is  a  function  of  particle  spatial  location,  particle 
velocity,  particle  mass,  particle  temperature  and  time.  The  trans¬ 
port  equation  for/,  which  is  derived  from  Boltzmann-BKG  approx¬ 
imation  model  of  gas  dynamics  [146],  is  given  by 

l  +  lt(fV)  +  l^A)=‘^rr  (76> 

where  A  is  the  particle  acceleration,  f  d  is  the  PDF  for  the  local  mass 
averaged  particle  velocity,  and  td  is  the  collision  damping  time. 
The  particle  acceleration  A  is  given  by 

,  dv  c  1  1 

A  =  —j—  =  9(vf-vs) - Vp+g - Vts  (77) 

dt  ps  asps 


where  vs  is  the  solid  velocity,  ps  is  the  solid  density,  as  is  the  solid 
volume  fraction,  ts  is  the  solid  normal  stress  and  9  is  the 
momentum  transfer  coefficient. 

The  equation  for  solid  movement  is  given  by 


(78) 


The  energy  equation  for  solid  is  given  by 


dTs  _  1  kfNUf  p 
v  dt  ~  ms  2 rs 


As(Tf  —  Ts) 


(79) 


where  C„  is  the  specific  heat  of  the  solid  material,  Ts  is  the  solid 
temperature,  Nufp  is  the  Nusselt  number  of  heat  transfer  in  the 
fluid  to  the  particle. 

The  particle  volume  fraction  as  is  defined  as 


«s 


dms  d~vs  dTs 


(80) 


jt(otfPfhf)+\'  ■  (afPfhf  Vf)  =  af^+Atf  ■  Vp'j  +®  +  Q+Sh  +  qD  +  V  •  ( afkfVTf ) 

(71) 

where  is  the  fluid  phase  enthalpy,  Tf  is  the  fluid  temperature,  0 
is  the  viscous  dissipation,  Q  is  the  energy  source  per  volume,  Sh  is 
the  energy  exchange  from  the  solid  face  to  the  fluid  phase,  and  qD 
is  the  enthalpy  diffusion.  The  thermal  conductivity  kf  is  the 
molecular  conductivity  plus  eddy  conductivity.  For  the  enthalpy, 
ideal  gas  approximation  is  used,  where  the  enthalpy  depends  on 
temperature  and  the  enthalpy  for  fluid  species  i  is  given  by 


From  conservation  of  volume,  the  sum  of  fluid  and  particle  volume 
fractions  equals  unity. 

af  +  as  =  1.  (81) 


The  fluid  momentum  equation  implicitly  couples  fluid  and  parti¬ 
cles  through  the  inter-phase  momentum  transfer.  The  inter-phase 
momentum  transfer  at  is  given  by 


F  = 


dms 
s  dt 


dm s  d~v  s  dTs 


(82) 


The  heat  transfer  from  the  solid  phase  to  the  gas  phase  is  given  by 


hi  = 


f  6p(  dt+Ahf  i 

JTrPf 


(72) 


where  Cp  is  the  specific  heat  at  constant  pressure  and  A hf  is  the 
heat  of  formation  of  the  species.  The  total  fluid  enthalpy  is  the 
mass  fraction  weighted  sum  of  the  fluid  species  enthalpies. 
Prandtl  number  correlation  is  given  by 


dms  dv  s  dTs 
(83) 


where  hs  is  the  particle  enthalpy. 


5.  Numerical  solution 


Prt  =  ^  (73) 

The  enthalpy  diffusion  term  is  given  by 

Qd  =  2  V  •  (hiafPfDYf  i)  (74) 

i  =  1 

The  mixture  properties  are  based  on  the  mass  fraction  of  each  gas. 
The  gas  phase  pressure,  temperature  and  density  are  related 
through  ideal  gas  equation.  For  an  ideal  gas,  the  pressure  is  given 
by  the  equation  of  state  as 

P  =  PfRTfl ^  (75) 


In  CPFD,  the  conservation  equations  are  integrated  over  a 
control  volume.  The  mixture  fluid  density,  velocity  and  pressure 
are  coupled  by  semi-implicit  pressure  equation  derived  from  the 
gas  mass  conservation  equation  which  is  applicable  to  an  arbitrary 
Mach  number.  The  fluid  momentum,  energy  and  pressure  equa¬ 
tions  are  solved  with  a  conjugated  gradient  solver.  The  chemistry 
of  differential  equations  are  calculated  using  a  stiff,  sparse  ODE 
solver. 

In  the  MP-PIC  scheme,  solids  properties  are  mapped  to  and 
from  the  Eulerian  grid  to  get  grid  properties  for  the  solids.  Fluid 
properties,  are  in  turn,  mapped  to  discrete  particle  locations.  The 
interpolation  operator  is  the  product  of  interpolation  operators  in 
the  three  orthogonal  directions.  For  a  particle  located  at  xs,  where 
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Fig.  13.  Gasifier  parameters  at  different  axial  levels  in  the  reactor  [149], 


xs  =  (xs,ys,zs),  the  x-directional  component  of  the  interpolation 
operator  to  grid  cell  i,  is  an  even  function,  independent  of  the  y 
and  z  coordinates,  and  has  the  properties. 


The  y  and  z  interpolation  operators  have  a  similar  form.  The 
particle  volume  fraction  at  cell  £  from  mapping  particle  volume  to 
the  grid  is 


Sj(xs) 


0,  x{_t  >xs>xf+I, 
1.  xs=xe 


<84>  =  nskSs(k 

k  Psfk 


(85) 
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Fig.  14.  (A)  Discrete  solids  colored  by  volume  fraction  mapped  from  the  grid  to  particle  locations  and  (B)  solids  field  colored  by  volume  fraction  [152], 


where  V4  is  the  grid  volume,  ns  is  the  number  of  particles  in  a 
numerical  particle  each  containing  a  cloud  of  particles  with 
identical  particle  velocity  us,  mass  ms,  temperature  Ts  located  at 
position  xs  and  the  summation  is  over  all  numerical  particles  Ns. 

The  inter-phase  drag  coefficient  is  given  by 


The  particle  normal  stress  is  an  approximation  of  collective  effects 
of  neighbor  particles  on  a  particle.  The  particle  stress  is  derived 
from  particle  volume  fraction  which  is  in  turn  calculated  from 
particle  volume  mapped  to  the  grid.  The  particle  normal  stress 
model  used  here  is  from  Harris  and  Crighton  [147]  as  given  below: 


e=cd 


3pf\vf-vs\ 
8  ft  ri 


(86) 


where  Cd  is  the  drag  coefficient.  The  correlation  for  Cd  based  on 
Wen  and  Yu  [77]  is  given  by 


JH(l  +  0.15Re0687); 
Q=\0.44, 

Re  < 1000 

Re  > 1000 

(87) 

and 

Pc  2pf\vf-v,\r, 

Pf 

(88) 

where  rs  is  the  particle  radius  as  given  below. 


rs  = 


m,a 


(89) 


s  max[(dCP-as),  f(l-as)] 

where  Ps  is  a  positive  constant  that  has  the  unit  of  pressure,  aCP  is 
the  solid  volume  fraction  at  close  packing  limit,  e  is  a  small 
number  on  the  order  of  10  - 7  to  remove  the  singularity  and  the 
constant  f  is  recommended  as  2  <  £  <  5. 

The  detailed  equations  and  solution  algorithm  for  Euler- 
Lagrange  CPFD  model  can  be  found  in  the  literature  [143], 

Snider  and  Banerjee  [148]  applied  CPFD  numerical  to  simulate 
the  ozone  decomposition  in  a  bubbling  fluidized  bed.  Simulations 
were  run  for  the  full  three-dimensional  bubbling  bed  of  0.229  m 
diameter  and  2  m  height.  The  ozone  decomposition  was  described 
by  a  single  stoichiometric  equation  with  first  order  reaction  rate 
and  isothermal  simulation  was  run.  The  simulated  ozone  mass 
fraction  as  function  of  inlet  velocity  was  compared  with  the 
analytical  solution  as  well  as  the  experimental  result  and  agreed 
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well  with  the  both.  Snider  et  al.  [149]  extended  the  Eulerian- 
Lagrangian  CPFD  methodology  to  include  the  heat  transfer  and 
chemistry  with  solid  material  pyrolysis.  The  model  was  used  to 
simulate  a  large  three-dimensional  fluidized  bed  coal  gasifier 
(4.57  m  diameter  and  13.7  m  height)  with  internal  cyclones, 
internal  heaters  and  as  sparger  to  illustrate  the  complexity  of 
problem  which  can  be  solved  by  the  CPFD  method.  The  gasifier 


was  initially  filled  with  47,000  kg  of  solids.  The  solid  was  modeled 
as  pure  carbon  with  a  particle  size  distribution.  The  gasification 
chemistry  was  described  with  five  homogeneous  and  heteroge¬ 
neous  chemical  reactions.  The  gas-solid  chemistry  consumed  the 
solids  and  solid  shrink  from  chemistry.  The  simulation  was  run  for 
300  s  which  took  6  days  computational  time  on  a  single  Intel® 
Core™  i7  computer.  The  results  from  the  simulation  provided 


Fig.  15.  Comparison  between  simulated  and  experimental  dry  gas  compositions:  (A)  T=800  C,  FR— 0.35,  S/B= 0.5;  (B)  T=850  C,  FR— 0.35,  F/B— 0.5;  (C)  T=800  C,  FR— 0.35, 
S/B  — 0.8;  (D)  7=800  C,  ER= 0.35,  S/B= 0.2;  (E)  7=800  °C,  ER=0.3,  S/B  =  0.5;  and  (F)  7=800  °C,  ER=0.4,  S/B= 0.5  [1521. 
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insight  into  the  gasifier  as  discussed  below.  Fig.  13  shows  that  the 
gas-solid  flow  had  regions  of  high  and  low  volume  fraction  and 
large  void  structures  (bubbles)  into  the  bed.  The  bed  surprisingly 
has  a  uniform  temperature  except  at  the  feed-sparger  arms  where 
high  sensible  heat  was  generated  from  combustion.  There  were 
variations  in  gas  species  radially  and  axially  in  the  bed  such  as 
pockets  of  H2  above  the  sparger-distribution-pipe.  Abbasi  et  al. 
[150]  simulated  the  feeding  section  of  a  fast  fluidized  bed  coal 
gasifier  by  using  CPFD  methodology.  The  solid  particles  are 
assumed  to  remain  unchanged  in  size  and  the  same  size  distribu¬ 
tion  was  taken  for  carbon  and  ash  particles.  Distribution  of  particle 
velocities,  particle  volume  fraction,  gas  compositions  and  tem¬ 
perature  were  obtained.  The  early  sing  of  chocking  in  the  feeding 
section  was  predicted  from  the  simulated  particle  volume  fraction 
distribution.  To  compare  the  CPFD  results,  a  simplified  approach, 
the  plug  flow  reactor  (PFR)  model  was  simulated  with  the  same 
operating  conditions.  It  was  found  that  the  time  and  cross-section 
average  value  of  temperature  and  gas  compositions  simulated  by 
CPFD  closely  approximated  the  one  for  PFR.  The  CPFD  model  was 
also  used  to  study  the  biomass  gasification  in  three-dimensional 
lab-scale  fluidized  bed  gasifier  [151,152]  by  considering  different 
particle  size  for  sand  and  biomass  particles.  The  gas-solid  flow 
pattern,  gas  composition  distribution,  parametric  variation  of  gas 
composition  etc.  were  investigated  in  details.  The  bubble  forma¬ 
tion,  development  and  eruption  with  time  was  well  captured  by 
complex  three  dimensional  flow  structures  revealed  by  the  simu¬ 
lation  as  shown  in  Fig.  14  [152]  whereas  only  the  growth  was 
observed  but  coalescence  and  eruption  was  not  identified  because 
of  high  velocity  [151],  It  was  also  observed  that  the  overall 
dynamics  of  the  bed  was  strongly  influenced  by  the  biomass 
injection  through  the  side  port  where  the  instantaneous  devola- 
tization  of  biomass  was  assumed.  The  mode  also  captured  the 
inherent  unsteady  nature  of  the  fluidized  bed.  The  comparison  of 
simulated  outlet  gas  composition  with  experimental  data  showed 
reasonable  good  agreement  (Fig.  15).  Nevertheless,  the  CPFD 
simulation  slightly  over-predicts  H2,  CO  and  CH4  generation  and 
under-predicts  C02.  It  was  concluded  that  the  little  discrepancies 
between  experimental  and  simulated  gas  compositions  may  be 
due  to  the  selection  of  kinetic  parameters  based  on  literature.  Xie 
et  al.  [153]  also  simulated  the  same  lab-scale  bubbling  fluidized 
bed  coal  gasifier  by  using  Euler-Lagrange  CPFD  model  which  was 
earlier  simulated  by  using  Euler-Euler  model  [89,91],  In  the 
present  model,  a  particle  size  distribution  was  considered  instead 
of  a  mean  particle  diameter.  The  model  considered  the  solid  mass 
consumed  or  produced  in  reactions  changed  the  size  of  particles. 
The  model  predicted  the  gas-solid  flow  pattern,  distribution  of  gas 
composition,  heterogeneous  and  homogeneous  chemical  reaction 


Fig.  16.  Cumulative  consumption  of  carbon  mass  versus  time  (Qc=8.0  l<g/h, 
Q^  — 28.4  kg/h,  <1=4.6  kg/h  [153]). 


rates  and  consumption  of  carbon  mass.  Fig.  16  shows  the  cumu¬ 
lative  net  consumption  of  carbon  mass  with  time.  It  was  described 
that  the  cumulative  consumption  of  carbon  particle  mass  over 
32.5  s  was  about  0.0098  kg/s,  which  resulted  in  the  average 
change  of  0.268  pm  for  solid  particle  radius.  The  simulated  outlet 
gas  composition  also  showed  good  agreement  with  the  experi¬ 
ment  with  minimum  relative  error  of  less  than  3%,  the  maximum 
relative  error  of  about  30%  and  the  average  relative  error  of  less 
than  15%. 

The  above  studies  showed  that  the  Euler-Lagrange  CPFD  model 
can  simulate  the  fluidized  bed  gasification  process  in  lab-scale  as 
well  as  commercial-scale  three  dimensional  gasifier  geometry 
with  reasonably  less  simulation  time  compared  to  the  Euler- 
Lagrange  DEM/DPM  model  and  also  by  keeping  the  Lagrangian 
description  of  particles  but  with  the  assumption  of  multiple 
particles  in  a  cell.  It  can  provide  many  details  inside  the  gasifier 
like  particle  shrinkage  due  to  reaction,  residence  time  of  particles, 
particle  size  distribution,  particle  attrition,  particle  agglomeration 
etc.  which  cannot  be  obtained  by  Euler-Euler  model.  It  also 
provides  the  three  dimensional  gas-solid  flow  structure,  gas 
composition  distribution,  temperature  distribution,  reaction  rate 
distribution  etc.  But,  the  Euler-Lagrange  CPFD  model  cannot 
provide  the  information  about  each  and  individual  particle  but 
can  provide  the  information  about  each  numerical  particle  (i.e. 
group  of  particles  having  similar  property)  and  hence,  it  is  a 
tradeoff  between  the  accuracy  and  the  computational  cost.  Differ¬ 
ent  models  used  to  simulate  the  fluidized  bed  gasification  process 
are  summarized  in  Table  8. 


6.  Conclusion 

An  updated  survey  of  published  mathematical  models  for 
fluidized  bed  gasification  process  is  presented.  The  present  paper 
describes  different  modeling  approaches  starting  from  the  simple 
models  like  equilibrium  model,  two-phase  flow  model  to  the  very 
complex  Euler-Euler  model  and  Euler-Lagrange  models.  The  basic 
principle  of  modeling  is  discussed  in  details.  Different  modeling 
approaches  used  by  different  researchers  are  reviewed  and  major 
results  obtained  by  them  are  also  discussed.  The  following  con¬ 
clusions  may  be  drawn  from  the  present  review. 

The  equilibrium  model  is  computationally  inexpensive.  It  does 
not  depend  on  the  gasifier  type.  But,  it  provides  a  quick  idea  about 
the  limiting  behavior  of  a  gasification  system  and  hence  it  is  useful 
for  preliminary  design  of  the  gasifier.  It  is  observed  that  the  results 
from  the  equilibrium  model  can  be  improved  by  modifying  it 
kinetically.  While,  in  kinetic  model  detailed  gasification  reaction 
kinetic  is  considered  with  respective  rate  of  reactions.  Kinetic 
models  also  consider  the  hydrodynamics  of  the  bed.  Depending  on 
the  hydrodynamic  consideration,  the  gasification  model  may  be 
classified  as  the  two-phase  flow  model,  the  Euler-Euler  model  and 
the  Euler-Lagrange  model.  The  two-phase  flow  modeling  is  the 
simplest  amongst  the  rate  based  or  kinetic  model  because  it  does 
not  solve  any  momentum  equation  for  gas  or  solid.  It  considers  the 
hydrodynamics  of  the  fluidized  bed  gasifier  which  is  based  on 
empirical  correlations.  In  Euler-Euler  model,  both  gas  and  solid 
are  treated  as  interpenetrating  continuum  and  separate  governing 
equations  are  solved  for  each  phase.  Coupling  is  achieved  through 
the  inter-phase  transfer  coefficients.  It  is  computationally  less 
expensive  compared  to  the  Euler-Lagrange  model  and  also  pro¬ 
vides  reasonably  good  details  about  the  fluidized  bed  gasifier. 
Euler-Lagrange  model  provides  the  maximum  number  of  details 
but  computationally  more  expensive  when  compute  individual 
particle  dynamics  (DEM/DPM).  Whereas,  the  Computational  Par¬ 
ticle  Fluid  Dynamics  (CPFD)  model,  is  a  modified  Euler-Lagrange 
model  which  assumes  a  numerical  particle  by  combining  particles 


Table  8 

Different  models  for  simulating  the  fluidized  bed  gasification  process. 


Authors 

Year 

Model  type 

Key  features 

Loha  et  al.  [6] 

2011 

Equilibrium  model 

Stoichiometric,  considered  only  gaseous  species,  modified 
equilibrium  constants. 

Loha  et  al.  [7] 

2011 

Equilibrium  model 

Stoichiometric,  considered  gaseous  and  solid  species. 

Li  et  al.  [10] 

2001 

Equilibrium  model 

Non-stoichiometric,  RAND  algorithm,  considered  gaseous 
and  solid  species,  modified  with  experimental  carbon 
conversion. 

Li  et  al.  [11] 

2004 

Equilibrium  model 

Non-stoichiometric,  RAND  algorithm,  considered  gaseous 
and  solid  species,  modified  with  experimental  carbon 
conversion  and  methane  formation. 

Petersen  and  Werther 
[22] 

2004 

Equilibrium  model 

Stoichiometric,  considered  only  gaseous  species. 

Srinivas  et  al.  [23] 

2009 

Equilibrium  model 

Stoichiometric,  considered  only  gaseous  species. 

Fryda  et  al.  [24] 

2008 

Equilibrium  model 

Non-stoichiometric,  fixed  amount  of  un-reacted  char  and 
CH4  and  tar  were  introduced  based  on  literature. 

Schuster  et  al.  [25] 

2001 

Equilibrium  model 

Non-stoichiometric,  considered  gaseous  and  solid  species. 

Petersen  and  Werther 
[22] 

2005 

Two-phase  flow 
model 

1.5-D,  unsteady-state 

Petersen  and  Werther 

[37] 

2005 

Two-phase  flow 
model 

3-D,  unsteady-state 

Jiang  and  Morey  [33] 

1992 

Two-phase  flow 
model 

1-D,  steady-state,  non-isothermal 

Chatterjee  et  al.  [40] 

1995 

Two-phase  flow 
model 

1-D,  steady-state,  non-isothermal 

Yan  et  al.  [41  ] 

1998 

Two-phase  flow 
model 

1-D,  steady-state,  isothermal,  net-flow  consideration 

Jennen  et  al.  [42] 

1999 

Two-phase  flow 
model 

1-D,  non-isothermal,  unsteady-state 

Hamel  and  Krumm 

[43] 

2001 

Two-phase  flow 
model 

1-D,  steady-state,  non-isothermal 

Fiaschi  and  Michelini 
[44] 

2001 

Two-phase  flow 
model 

1-D,  non-isothermal 

Sadaka  et  al.  [45-47] 

2002 

Two-phase  flow 
model 

1-D,  unsteady-state,  non-isothermal 

Chejne  and 

Hernandez  [48] 

2002 

Two-phase  flow 
model 

1-D,  steady-state,  non-isothermal. 

Ross  et  al.  [49] 

2005 

Two-phase  flow 
model 

1-D,  steady-state,  non-isothermal. 

Radmanesh  et  al.  [50] 

2006 

Two-phase  flow 
model 

1-D,  steady-state,  Isothermal 

Kaushal  et  al.  [51] 

2010 

Two-phase  flow 
model 

1-D,  steady-state,  non-isothermal 

Application 


Output 


BFB,  lab-scale,  biomass  gasification,  steam 
as  gasifying  agent. 

BFB,  lab-scale,  biomass  gasification,  steam 
as  gasifying  agent. 

CFB,  lab-scale,  coal  gasification,  air  as 
gasifying  agent. 


Parametric  variation  of  gas  composition,  correlation  for  product  gas. 

First  law  and  second  law  efficiencies,  HHV,  external  energy  input, 
CBP. 

Parametric  variation  of  gas  composition,  HHV,  carbon  conversion, 
carbon  formation  boundary,  role  of  moisture 


CFB,  lab-scale,  coal  gasification,  air  as 
gasifying  agent. 


Parametric  variation  of  gas  composition,  gas  yield,  HHV,  cold  gas 
efficiency,  carbon  formation  boundary. 


CFB,  pilot-scale,  sewage  sludge 
gasification,  air  as  gasifying  agent. 

CFB,  pressurized  condition,  biomass 
gasification,  air-steam  as  gasifying  agent. 
FB,  lab-scale,  biomass  gasification,  steam 
as  gasifying  agent. 

Dual  fluidized  bed,  commercial-scale, 
biomass  gasification,  steam  as  gasifying 
agent. 


Parametric  variation  of  gas  composition. 

Parametric  variation  of  gas  composition,  heating  value,  exergy, 
efficiency. 

Parametric  variation  of  gas  composition,  efficiency,  integrated  with 
SOFC 

Parametric  variation  of  gas  composition,  heating  value,  efficiency, 
carbon  formation  boundary. 


CFB,  pilot-scale,  sewage  sludge 
gasification,  air  as  gasifying  agent 
CFB,  pilot-scale,  sewage  sludge 
gasification,  air  as  gasifying  agent. 

BFB,  lab-scale,  biomass  gasification,  air  as 
gasifying  agent. 

BFB,  lab-scale,  coal  gasification,  steam  and 
air  as  gasifying  agent. 

BFB,  coal  gasification. 


CFB,  pilot-plant,  biomass  gasification,  air 
as  gasifying  agent. 

BFB,  lab-scale  to  commercial-scale. 
Gasification  of  brown  coal,  peat  and 
sawdust.  Air,  air/steam  and  02/steam  as 
gasifying  agent. 

BFB,  biomass  gasification,  air  as  gasifying 
agent. 

Dual  distributor  type  FB,  lab-scale, 
Biomass  gasification.  Air-steam  as 
gasifying  agent. 

BFB,  lab-scale  coal  gasification. 


BFB,  commercial-scale  and  pilot-scale, 
coal  gasification,  air-steam  as  gasifying 
agent. 

BFB,  biomass  gasification,  air-steam  as 

gasifying  agent 

BFB,  biomass  gasification 


Axial  gas  composition  variation,  comparison  with  own  experimental 
data,  modified  the  reaction  kinetics. 

Three  dimensional  variation  of  product  gas,  uneven  distribution 
near  feeding  point,  better  mixing  with  by  increasing  number  of 
feeding  point. 

Fuel  feed  rate,  gas  composition,  heating  value,  temperature  etc. 
model  agreed  well  with  experiment  at  higher  temperature  but  failed 
at  lower  temperature. 

Variation  of  gas  composition,  temperature,  calorific  value  and 
carbon  conversion  with  oxygen  feed  and  steam  feed. 

Bed  voidage,  reaction  rate  and  gas  composition  along  the  height, 
carbon  conversion,  significant  deviation  without  considering  the 
net-flow. 

Axial  gas  composition  and  temperature  profile. 

Overall  carbon  conversion,  freeboard  temperature,  gas  composition. 


Gas  composition  and  temperature  variation  along  the  axis, 
optimization  with  respect  to  ER,  pressure,  bed  height  and  gas 
velocity. 

Gas  composition,  bed  temperature,  heating  value,  gas  production 
rate. 

Temperature,  gas  composition,  volume  fraction,  velocity  and  other 
fluid  dynamic  parameters.  Strong  influence  of  feed  point  position 
was  observed  in  the  result. 

Gas  composition,  temperature  and  reaction  rate  along  the  height, 
overall  carbon  conversion. 

Gas  composition  variation  with  T,  ER,  S/B,  feed  location  and  mass 
transfer  between  counter  current  region. 

Temperature,  solid  hold  ups  and  gas  concentration  variation  along 
the  reactor’s  major  axis. 
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Table  8  ( continued ) 


Authors 

Year 

Model  type 

Key  features 

Application 

Output 

Pengmei  et  al.  [52] 

2008 

Two-phase  flow 
model 

1-D,  steady-state,  Isothermal 

FB,  biomass  gasification,  air-steam  as 
gasifying  agent 

Gas  composition,  gas  yield 

Goyal  et  al.  [53] 

2010 

Two-phase  flow 
model 

1-D,  non-isothermal 

BFB,  gasification  of  coal  and  petcoke 
mixture 

Effect  of  composition  and  location  of  feed  point  and  ash  content  on 
the  performance,  Increase  in  petcoke  content  tends  to  lower  the 
efficiency  and  carbon  conversion  but  increases  the  amount  of  syngas 
produced,  increase  in  ash  content  of  coal  decreases  the  carbon 
conversion. 

Gungor  [54] 

2011 

Two-phase  flow 
model 

1-D,  steady-state,  isothermal 

FB,  biomass  gasification,  air-steam  as 
gasifying  agent. 

Hydrogen  production  variation  with  T,  ER,  S/B,  velocity 
and  particle  size. 

Dou  et  al.  [84] 

2008 

Euler-Euler  model 

2-D,  single  solid  phase  and  single  particle  size,  unsteady- 
state,  non-isothermal. 

BFB,  lab-scale,  steam  reforming  of 
glycerol. 

Gas-solid  flow  pattern,  glycerol  consumption,  hydrogen 
and  other  gas  production. 

Dou  and  Song  [85] 

2010 

Euler-Euler  model 

2-D,  single  solid  phase  and  single  particle  size,  unsteady- 
state,  non-isothermal. 

BFB,  lab-scale,  steam  reforming  of 
glycerol. 

Gas-solid  flow  pattern,  glycerol  consumption,  hydrogen 
and  other  gas  production,  detailed  analysis  of  relationship 
between  hydrodynamics  and  hydrogen  production. 

Papadikis  et  al. 

[86,87] 

2008 
&  2009 

Euler-Euler  model 

2-D  and  3-D,  Eulerian  description  of  gas  and  sand  particle, 
Lagrangian  description  of  single  biomass  particle. 

BFB,  Lab-scale,  Fast  pyrolysis 

Gas-solid  flow  pattern  with  biomass  particle  position,  gas 
composition  and  temperature  distribution,  tar  evolution 

Xue  et  al.  [88] 

2012 

Euler-Euler  model 

2-D  and  3-D  simulation,  two  solid  phases  for  sand  and 
biomass,  unsteady-state,  isothermal. 

BFB,  lab-scale,  biomass  fast  pyrolysis 

Influence  of  operating  conditions  on  the  conversion  process  and 
product  yield,  model  over-predicted  the  un-reacted  biomass 
elutriation  for  smaller  particles. 

Yu  et  al.  [89] 

2007 

Euler-Euler  model 

2-D,  k-e  turbulence  model  for  gas  phase,  Single  solid  and 
single  particle  size,  unsteady-state,  non-isothermal. 

BFB,  lab-scale,  coal  gasification,  air-steam 
as  gasifying  agent. 

Gas-solid  flow  pattern,  gas  composition  distribution,  reaction  rate 
distribution. 

Wang  et  al.  [90] 

2009 

Euler-Euler  model 

3-D,  k-e  turbulence  model  for  gas  phase,  Single  solid  phase 
and  single  particle  size,  unsteady-state,  non-isothermal. 

BFB,  lab-scale,  coal  gasification,  air-steam 
as  gasifying  agent. 

Flow  pattern,  gas  velocity,  particle  velocity,  gas  composition 
distribution,  reaction  rate  distribution. 

Armstrong  et  al.  [91] 

2011 

Euler-Euler  model 

2-D,  separate  solid  phase  for  coal,  limestone  and  char.  Single 
particle  size  for  each  solid  phase,  unsteady-state,  non- 
isothermal. 

BFB,  lab-scale,  coal  gasification,  air-steam 
as  gasifying  agent. 

Gas-solid  flow  pattern,  distribution  of  gas  composition,  reaction  rate 
and  temperature.  Impact  of  limestone  calcinations  on  gaseous 
composition. 

Armstrong  et  al.  [92] 

2011 

Euler-Euler  model 

2-D,  two  different  devolatilization  model,  separate  solid 
phase  for  coal,  limestone  and  char,  single  particle  size  for 
each  solid  phase,  unsteady-state,  non-isothermal.  extended 
simulation  time  (400  s). 

BFB,  lab-scale,  coal  gasification,  air-steam 
as  gasifying  agent. 

Gas-solid  flow  pattern,  gas  composition  variation  with  bed  material 
composition,  bed  height,  temperature  of  bed  and  inclusion  of  heat 
transfer  coefficient. 

Silaen  and  Wang 

[138] 

2010 

Euler-Lagrange 

DEM/DPM 

3-D,  four  different  devolatization  models  are  employed 

Entrained  flow  gasifier,  coal  gasification, 
oxygen  blown. 

Gas  composition,  temperature,  efficiency,  effect  of  different 
devolatization  model. 

Watanabe  and  Otaka 

[139| 

2006 

Euler-Lagrange 

DEM/DPM 

2-D,  2  t/day 

Entrained  flow  gasifier,  coal  gasification. 

Influence  of  the  air  ratio  on  carbon  conversion  efficiency,  amount  of 
product  char,  heating  value,  and  cold  gas  efficiency,  temperature  and 
product  gas  composition  distribution. 

Grabner  et  al.  [140] 

2007 

Euler-Lagrange 

DEM/DPM 

4800  t/day,  high  pressure  (33  bar). 

Winkler  gasifier,  coal  gasification,  steam/ 
oxygen  as  gasifying  agent. 

Flow  pattern,  turbulence,  product  gas  composition,  temperature  and 
radial  heat  transfer. 

Bruchmuller  et  al 

[141] 

2012 

Euler-Lagrange 

DEM/DPM 

3-D,  0.8  million  individual  sand  and  biomass  particle. 

BFB,  lab-scale,  biomass  fast  pyrolysis. 

Gas-solid  flow  pattern,  particle  degradation,  particle  entrainment, 
gas  composition  distribution. 

Snider  and  Banerjee 
[148] 

2010 

Euler-Lagrange 

CPFD 

3-D,  constant  particle  size,  only  single  reaction,  unsteady- 
state. 

BFB,  lab-scale,  ozone  decomposition. 

Particle  volume  fraction  distribution,  ozone  mass  fraction 
distribution,  variation  with  velocity. 

Snider  et  al.  [149] 

2011 

Euler-Lagrange 

CPFD 

3-D,  47,000  kg  of  particle  initially,  particle  size  distribution, 
particle  shrinkage  from  chemistry,  unsteady-state,  non- 
isothermal. 

Internal  CFB,  commercial-scale  example 
gasifier,  complex  geometry  with  internal 
cyclones,  coal  gasification,  air-steam  as 
gasifying  agent. 

Particle  volume  fraction,  gas  composition,  temperature  and  reaction 
rate  distribution. 

Abbasi  et  al.  [150] 

2011 

Euler-Lagrange 

CPFD 

2-D,  same  particle  size  distribution  for  coal  and  ash  particle, 
size  does  not  change  due  to  reaction,  unsteady-state,  non- 
isothermal. 

Feeding  section  of  a  CFB,  lab-scale,  coal 
gasification,  air-steam  as  gasifying  agent. 

Particle  volume  fraction,  and  size  distribution,  fluid  velocity, 
temperature  and  pressure  drop  distribution,  gas  composition 
distribution. 

Xie  et  al.  [151  ] 

2012 

Euler-Lagrange 

CPFD 

3-D,  different  particle  species  for  carbon,  ash  and  sand 
particle,  unsteady-state,  isothermal. 

BFB,  lab-scale,  biomass  gasification, 
air-steam  as  gasifying  agent. 

Bubble  growth  was  captured  but  collision  and  eruption  were  not 
visible,  gas  composition  variation  with  T,  S/B,  ER.  Gas  composition 
distribution. 

Loha  et  al.  [152] 

2014 

Euler-Lagrange 

CPFD 

3-D,  different  particle  species  for  carbon,  ash  and  sand 
particle,  particle  size  distribution,  unsteady-state,  isothermal. 

BFB,  lab-scale,  biomass  gasification,  steam 
as  gasifying  agent. 

Bubble  formation,  development  and  eruption  with  time  and 
interaction  with,  influence  of  biomass  injection  through  side  port, 
gas  composition  variation  with  T,  S/B,  ER  and  time. 

Xie  et  al.  [153] 

2013 

Euler-Lagrange 

CPFD 

3-D,  different  particle  species  for  carbon,  ash  and  sand 
particle,  particle  size  distribution,  particle  size  changes  due 
to  reaction,  unsteady-state,  non-isothermal. 

BFB,  lab-scale,  coal  gasification,  air-steam 
as  gasifying  agent. 

Bubble  formation,  development  and  eruption  with  time,  distribution 
of  gas  composition,  temperature  and  reaction  rate,  carbon 
consumption  with  time. 
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of  same  properties,  can  simulate  large-scale  fluidized  bed  gasifier 
with  comparatively  less  time.  But,  the  information  about  each 
individual  particle  cannot  be  provided  by  the  Euler-Lagrange 
CPFD  model. 

From  the  present  survey,  it  is  observed  that  almost  in  all 
fluidized  bed  gasification  models  only  the  outlet  gas  compositions 
are  compared  with  experiment.  There  are  few  measurements 
available  for  comparison  with  detailed  model  results.  Therefore, 
effort  is  still  required  to  validate  the  truly  comprehensive  fluidized 
bed  models. 
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